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ABSTRACT 


This  thesis  investigates  the  viability  of  a  new  method 
for  numerically  computing  the  input  impedance  and  the  currents 
on  simple  antenna  structures.  This  technique  considers  the 
antenna  between  two  ground  planes  and  uses  multiregion 
cylindrical  harmonic  expansions  with  tangential  field  con¬ 
tinuity  to  obtain  the  surface  currents  and  input  impedance. 
The  computed  results  are  compared  to  the  results  obtained  from 
the  Numerical  Electromagnetics  Code  for  various  physical 
parameters  to  assess  computational  accuracy. 
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I.  INTRODUCTION 


There  are  a  number  of  methods  in  use  today  that  numeri¬ 
cally  solve  for  the  currents  on  simple  wire  antenna  struc¬ 
tures.  The  majority  of  these  methods  utilize  an  integral 
equation  solution  to  Maxwell’s  Equations  which  also  match  the 
boundary  conditions  on  the  structure  and  satisfy  the  radiation 
conditions  for  the  radiated  fields.  One  of  the  most  widely 
used  implementations  of  the  integral  equation  method  is  the 
Numerical  Electromagnetics  Code  (NEC)  [Ref  6).  This  program 
has  been  proven  to  be  highly  accurate  for  antenna  structures 
that  can  be  modeled  using  individual  wire  segments  or  surface 
patch  elements. 

A  new  method  is  presented  here  which  does  not  rely  upon 
an  integral  equation  formulation.  This  technique  encloses  the 
antenna  structure  between  two  parallel  ground  planes,  as  an 
approximation  to  the  free  space  condition,  to  allow  a  partial 
field  representation  using  a  periodic  Fourier  series.  In  the 
initial  investigation  being  reported  here,  the  simple  monopole 
radiator  is  being  considered.  In  such  a  case,  cylindrical 
regions  are  established  around  the  structure  and  the  harmonic 
expansions  of  the  magnetic  and  electric  fields  are  matched 
along  the  regional  boundaries  to  solve  for  the  unknown 
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expansion  coefficients.  These  coefficients  are  then  used  to 
compute  the  current  imposed  on  the  surface  of  the  antenna. 

This  project  began  as  a  search  for  the  viability  of 

computing  the  currents  on  a  top  loaded  monopole  between  two 
ground  planes  using  cylindrical  harmonic  expansions  in  three 
regions.  After  a  number  of  weeks  of  inconclusive  results,  it 
was  decided  to  look  at  the  simpler  case  of  the  monopole 

antenna  between  two  ground  planes  as  presented  in  Chapter  2. 

During  the  validation  of  the  code,  as  presented  in 

Chapter  3,  it  was  found  that  the  number  of  terms  used  in  each 
region  can  be  reduced  by  as  much  as  90%  and  still  yield 
consistent  results  for  the  surface  current  and  input  im¬ 
pedance.  However,  the  representation  of  the  electric  field 
along  the  surface  of  the  antenna  is  highly  inaccurate  when  the 
number  of  terms  are  insufficient  to  adequately  represent  the 
field.  This  phenomenon  will  be  discussed  in  Chapter  4. 
Additional  research  revealed  that  a  simpler  method  could  be 
used  that  would  remove  the  requirement  for  one  complete 

expansion.  This  is  presented  in  Chapter  5.  Finally,  con¬ 
clusions  and  recommendations  are  presented  in  the  last 
chapter.  In  addition,  Appendices  are  used  to  contain  much  of 
the  detailed  descriptions  of  the  computer  algorithms  and 
resultant  programming  for  this  effort. 
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A.  METHOD  USED 


The  general  approach  is  to  place  the  antenna  structure 
being  analyzed  between  two  perfectly  conducting  ground  planes 
and  divide  the  surrounding  area  into  a  number  of  cylindrical 
regions,  using  a  cylindrical  coordinate  system,  that  coincide 
with  the  natural  boundary  interfaces  of  the  antenna  as  shown 
in  Figure  l.  The  upper  ground  plane  is  introduced  to  com¬ 
pletely  close  the  regions  and  produces  multiple  images,  each 
identical  to  the  previous  image  but  with  opposite  polarity. 

Once  the  regions  are  defined,  we  start  by  obtaining  the 
solution  to  the  complex  scalar  Helmholtz  equation.  Then,  by 


using  an  harmonic  expansion  of  the  resulting  scalar  wave 
function,  ifi,  we  write  an  expression  for  the  vector  electric 
and  magnetic  fields  in  each  region.  The  tangential  component 
of  the  electric  field,  as  approximated  by  the  truncated 
expansion  in  Region  I,  is  set  equal  to  the  "known”  field  in 
the  gap  and  set  equal  to  zero  along  the  cylindrical  surface 
of  the  perfectly  conducting  monopole.  The  tangential  electric 
and  magnetic  field  expansions  are  then  matched  across  the 
boundary  interface  between  Regions  I  and  II.  The  resulting 
set  of  equations  is  reduced  to  a  system  of  one  equation  with 
one  unknown  by  sifting  out  some  of  the  coefficients  using  the 
orthogonality  principle  of  Fourier  moment  integrations.  The 
resulting  system  is  then  solved  for  the  remaining  set  of 
unknown  coefficients.  These  coefficients  are  then  used  to 
compute  the  current  distribution  along  the  surface  of  the 
antenna.  To  validate  the  source  '■'ode,  the  results  are  then 
compared  to  those  obtained  from  the  Numerical  Electromagnetics 
Code  (NEC)  to  ascertain  computational  accuracy. 

B.  ASSUMPTIONS  MADE 

In  order  to  simplify  the  derivations,  assumptions  were 
made  regarding  the  electric  characteristics  along  the  ground 
planes  and  the  surface  of  the  antenna  structure,  and  the 
electric  field  produced  by  a  driving  source  of  constant 
frequency.  First,  the  ground  planes  are  modeled  as  perfectly 
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conducting  ground  planes  of  infinite  dimension.  In  reality, 
a  ground  plane  can  be  considered  to  have  infinite  dimension, 
for  purposes  of  antenna  current  calculations,  if  it  is  larger 
in  dimension  than  about  10  times  the  size  of  the  antenna  [Ref. 
1].  Additionally,  the  electric  field  produced  by  the  driving 
source  is  considered  to  be  of  constant  magnitude  within  the 
gap  region.  This  is  a  standard  source  model  employed  in 
integral  equation  modeling  of  antenna  structures,  [Ref.  1]. 
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II.  COMPUTATIONAL  PROCEDURE 


This  chapter  will  present  the  mathematical  background 
required  to  solve  for  the  currents  and  the  input  impedance  on 
a  simple  monopole  antenna  structure  using  cylindrical  harmonic 
expansions.  The  final  solution  will  be  reduced  to  a  system 
of  one  set  of  unknowns  which  will  be  solved  using  a  digital 
computer  system.  An  overview  of  the  software  code  used, 
including  error  detection,  hardware  requirements  and  a  brief 
discussion  of  the  input  and  output  file  structures  will  also 
be  presented.  Details  of  the  algorithm  implementation  are 
relegated  to  Appendices. 

A.  THEORETICAL  TECHNIQUE 

The  antenna  structure  shown  in  Figure  1  (page  3)  is 
orientated  within  the  cylindrical  coordinate  system  as  shown 
in  Figure  2.  We  begin  this  development  with  the  complex 
scalar  wave  equation,  known  as  the  Helmholtz  equation,  given 
in  cylindrical  coordinates  as  [Refs.  2,  3] 

2  2 
V  y/(p,z)  +  k  y(p,z )  =  0 

(Eq.  1-1) 

where  0  is  used  to  denote  the  electric  type  Hertz  potential 
and  not  the  magnetic  flux  [Ref.  3],  Expanding  this  equation, 
and  noting  that  the  structure  and  the  fields  are  assumed  to 
be  axisymmetric,  one  obtains  the  following: 
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Figure  2  Cylindrical  Coordinate  System 


l—[p 

P  dp  V  dP, 


2 

+i k  V  =0 

dz2 


{Eq.  1-2} 

Utilizing  axisymmetry,  the  transverse  magnetic  (TM) 
solutions  are  given  in  cylindrical  coordinates  as  [Ref.  3] 


E  «IiLsi 

p  y  dp  dz 


Hp  =0 


3y/ 


£z 


i  a 


i  0Z 


Y' 


where 


(Eq.  1-3} 


y  =jO)E0 


{Eq.  1-4} 
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The  following  boundary  conditionons  must  be  satisfied: 

•  The  electric  field  tangential  to  the  ground  planes  must 
be  zero  (Ep  =  0) . 

•  The  electric  field  tangential  to  the  top  of  the  antenna 
structure  must  be  zero. 

•  The  tangential  electric  field  along  the  antenna  structure 
must  be  zero  (Ec  =  0) . 

•  The  electric  field  along  the  gap  distance  d  must  be  a 
constant  value  (Ec  =  -V^d) . 

»  The  electric  and  magnetic  fields  must  be  continuous 
across  the  boundary  between  the  two  regions  above  the 
antenna  structures  (E,,  H*) . 

•  The  radiated  fields  in  Region  I  must  be  purely  outbound 
at  large  distances  from  the  antenna  and  must  satisfy  the 
radiation  condition  that  Et  =  r}0H*,  where  rjo  =  3770. 

Using  separation  of  variables,  where  the  product  solution 
has  the  form  =  R(p)Z(z),  one  obtains  the  following  equa¬ 
tions:  [Ref.  3] 


^+kt2Z  =0 

dz  (Eq.  1-5} 
The  first  equation  is  recognized  as  Bessel's  equation  of  order 
n.  The  solutions  to  Bessel's  equation  are  represented  by  any 
of  several  special  functions,  including  [Ref.  3] 


Jn (kpP)  ,YJkpP)  (Eg.  1-6} 
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where  J n(kpp)  and  Y a{kpp)  represent  Bessel  functions  of  the 

first  and  second  kind  of  order  n,  and  H na>(*p?)  and  H n<2>  (*„/>) 

are  Hankel  functions  of  the  first  and  second  kind  of  order  n, 
(2) 

where  H„  =  Jn  ±  j  Yn.  As  any  two  of  these  functions  are 
linearly  independent  of  each  other,  the  solution  to  Bessel's 
equation  can  be  represented  as  a  linear  combination  of  any  two 
of  these  functions.  The  solution  to  the  Z  equation  is  a 
linear  combinations  of  harmonic  functions  of  the  form  eJk*  and 
e'Jk* . 

Using  these  results,  the  desired  solutions  to  the  Helmho¬ 
ltz  equation  becomes 


where 


f 

P) 


(  sin  (kz  z)  ) 
jcos  {kzz)  J 


{Eq.  1-7} 


k 


2 

P 


2 

=  CO  HE 


(Eq.  1-8} 

and  where  k0  is  defined  as  the  wave  number,  co  the  angular 
frequency,  n  the  permeability  of  a  vacuum  (4jt  x  10';  H/m)  and 
c  the  permittivity  of  a  vacuum  (8.854  x  10’12  F/m) .  Acceptable 
quantized  values  of  kt  are  found  by  substituting  the  product 
solution  into  Equation  1-3  and  setting  E,  =  0  on  the  upper 
ground  planes  in  Region  I  while  enforcing  E,  =  0  on  the  upper 
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ground  plane  and  the  top  of  the  antenna  in  Region  II.  This 
gives,  respectively. 


K  - 

n 

nr 

/ 

in  Region  I 

(Eq.  1-9. a} 

K  = 

nr 

in  Region  II 

{Eq.  1-9. b} 

n 

IT 

It  should  be  noted  here 

that 

the  trigonometric 

terms  used  in 

Equation  1-7  will  uniquely  satisfy  the  first  three  boundary 
conditions  listed  previously  for  each  term  in  the  series. 

For  large  arguments,  the  Bessel  and  Hankel  functions 
listed  above  closely  resemble  harmonic  functions.  As  seen  in 
Table  I  [Ref.  3]/  J„ (kpp)  and  Ya(kpp)  are  analogous  to  cosine 
and  sine  functions  respectively,  while  Hn<1>(*pp)  and  H nlZ)  (kpp) 
are  related  to  complex  exponential  functions.  Therefore, 
Bessel  functions  can  be  used  to  represent  standing  waves, 
while  the  Hankel  functions  describe  travelling  waves;  specifi¬ 
cally  Hn<1>  represents  inbound  waves  while  H„(2)  represents 
outbound  waves.  Using  these  properties,  the  potential  field 
in  Region  I  can  be  written  as  the  product  of  Hankel  and  cosine 
functions  where  the  cosine  term  implicitly  satisfies  the 
tangential  electric  and  magnetic  field  boundary  conditions. 

Because  the  higher  spectral  bandwidth  of  the  potential 
field  generated  by  the  gap  voltage  was  expected  to  require  a 
larger  number  of  terms  in  the  Fourier  series  expansion  than 
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functions  •«  I  snd  Ioc  S.  »  -  0-  “d  '*  “d  *' 


that  of  the  quasi-singular  field  at  the  erd  of  the  antenna, 
the  total  field  in  Region  I  was  separated  into  two  parts  as 
shown  below 


w  wo 

Y(P>z)  =  X  o\vnp) cos  +  X  bmHo\vmp)cos 


n  =  0 

=  V a  +  n 


m  =  0 


/mttz 


/  ) 


{Eg. 1-10} 


where 


{Eg.  1-11} 

In  Region  II,  the  potential  field  is  constrained  as  a 
standing  wave  and  can  be  expressed  as 


where 


Y'CP.O*  X  cnJ0(unp)  cos 


n  *  0 


«;r 

lt 


M) 


{Eg.  1-12} 


{Eq.  1-13} 

and  q  =  l  -  h. 

Using  Equation  1-3,  the  electric  field  component  at  the  end 
of  the  antenna  is  approximated  by  the  truncated  expansion 
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Ni 


1  x. '  «  (2) 

Ez  (a,z)  =  x  >  vn  an  H 0V  '(vna)cos 

U  V 


^  n  =  0 


n  nz 

~r 


0, 


for 


N 2  1 

1  V*  2  r  /  \  n  ft(z  -  h) 

~  £iUncn  */0  (“„<*)  COS  ^ - - - j, 

y  n=  0 


for  h  <>  z  <,  l 


(Eq.  1-14} 

and  the  electric  field  due  to  the  gap  voltage  is  given  by 


M  i 

E zb {a,2  )  *  —  ^ 

y  m  =  0 


v„ 


„  (2).  v  [rmtz\ 
Ho  ^ma)cos\—T-J 


{-  ~  ,  for  0£z  £d 
0 ,  for  d  £  z  £  l 


{Eq.  1-15} 


where  Nj  represents  the  number  of  terms  to  be  used  to  in 
Region  I,  N2  for  Region  II,  and  Mx  for  the  electric  field 
component  due  to  the  gap  voltage. 

The  last  boundary  condition  to  enforce  is  that  the 
magnetic  field  be  continuous  across  the  interface  of  Region 
I  and  Region  II,  Again  using  Equations  1-8,  1-10  and  1-12, 
and  setting  H4J(a,z)  *  H^u(a,z),  an  expression  for  the  phi 
component  of  the  total  magnetic  field  can  be  written  as 
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", 

Hq  (a,z ) »  vn  an  H  }2)<yj i)  cos 

n  =o 


+  X  Vmb'»Hl2)lV»Py>C0S 


Ni 

^^uncnJl(urla)cos 


nx(z-h) 
<1  . 


n=0  {Eq.  1-16} 
Equation  1-16  can  be  reduced  to  one  set  of  unknowns  by  solving 
for  the  coefficients  ba  analytically  and  using  this  result  to 
solve  for  the  c„  in  terms  of  bn  and  the  unknown  t^'s.  This  is 
accomplished  by  sifting  out  the  desired  coefficient  using  the 
orthogonal  property  of  Fourier  moment  integrations  which  can 
be  defined  over  the  regions  where  the  fields  exist  as 


I 


(2) 

m 


m  *n 
m  =  n  * 0 
m  -  n  =0 

{Eq.  1-17. a} 

0,  m  *  n 

j,  m  =  n  *0 

q,  m  =n  =0 


{Eq.  1-17. b} 
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f  (m  n  z\ 

=  I  cosi— J 

J  n 


dz  = 


l  .  (m  K  d\ 

—  sm^  j  J, 


m  =  0 


m  *  0 


{Eq.  1-17. c} 


Applying  the  integrations  to  both  sides  of  Equation  1-15,  the 
bm  coefficients  can  be  solved  for  in  terms  of  known  functions 
resulting  in  an  expression  for  the  terms  in  the  second 
summation  of  Equation  1-16,  as  shown  below: 


M|  2  C  1 

2  bm  (Vma  )  I  COS  COS  l~ 

m-0  J  q  '  ’  ' 


-  -  y 


HKZ  \  , 
COS  — j — jdz 


{Eq.  1-18} 


Therefore, 


(Eq.  1-19) 


which  can  be  rearranged  to  give 


v-  «o<2,(v„0)  (Eq.  1.20j 

By  applying  the  Fourier  moment  integrations  to  Equation 
1-16  then  substituting  in  the  expression  for  ba  (Equation 
1-20) ,  the  cn  term  can  be  written  as 
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2  _  ,  .  J  o(una) 

M«  cnJ0(una)  =  un- 


w. 


^  1  (Mnfl)  I  *  =  0 


X  vkakHi2\vka)Tkin 


where 


J0£Yi  !hllT 

J  °  d  2a  y  ,  (1)  u  (2),  *•" 


— 0  v*  ir  H o  (v*  0 ) 


{Eq.  1-21) 


f1  (knz)  \nn  (  . 

*./»“  I  COS^-y— JC0S|~(Z  -/l)J 


dz 


{Eq.  1-22) 

We  now  have  one  equation  with  two  unknowns.  A  second  expres¬ 
sion  for  the  unknown  an's  can  be  obtained  by  applying  the 
moment  integrations  to  the  Equation  1-14  such  that 


C  1  i  \ 

r  \ 

1  cos  — 

cos 

m/rz 

1  cosi  /  / 
'  0 

i  / 1 

Nt  r  i 

=  X  UncnJ0^na)  I  COS 

«  =  0  J . 


\nK[z-h\\  lmnz\  . 

I — 5 - -  cos  - — ; — I  dz 

l  9  /  /  / 


{Eq.  1-23) 
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which  gives 


v„a„Ho 


Ni 

un  cit^O  (una)  T m  n 

n  =  0 


Rewriting  Equation  1-21  yields 


{Eq.1-24} 


N, 


UnCnJl(“na)In2)ss  V* °k ^  Tk.' 


k  =  0 


M, 

vo  V  1 

-'"Sr  2.  - 


Vn  ^  1  Ikm  Hjm(vka) 


(l)  u  <2), 


~ d  v*  Ik  H0"'(vka) 


k,n 


{Eq.  1-25) 


Substituting  this  result  into  Equation  1-24  yields 


N 


N 


»  02,(V„  «)I.0,-£  u„^;a)r^-  <  I  V,  a,  H  ,%  <■ )  r  M 

n  =0  'it  1  1  =0 


Vo  v  U(«na)TmJ  v  1  «12)(^ 

d  ZjUn  .(2).  ,  .  \  „  r  0)  „  (2),.. 


J°)E0  d  S-J  ""  .  (2).  ,  . 


v.  /  ^  //  <*Vv  d  1 
Jk  «0  *  /k  H  o  {Vka) 


k,n 


(Eq.  1-26) 


Collecting  terms  and  rewriting  yields 
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Iv^sw 


1=0  I'JiM 


-  vmfln// o  (vMa)/„ 


„yw  v°  y‘  [^±^1  y^Jo(una)TmMTkM 
0  d  k  =  0  VkhilW2)(Vk«)  1  „=0  /„(2)^  («»«) 


(Eq.  1-27} 


This  can  be  written  in  matrix  form  as 


where 


T,  Am,kak  =  Bm  for  m  =  0,  ! 


{Eq.  1-28} 


Am.k  =  vkH  \  \vka)Vm<k  m  *k 


Am,m  ~  vm  ^  1  (Vm  0  )  '  Vm  ^  )(vm  d 


m  '  m 


m  1  0  d  2-*  .  (i)  (2). 


k-oVklkHoWk*) 


{Eq.  1-29} 


,,-x 


unJo(una)TmnTkn 


n=  0 


{Eq.  1-30} 


The  current  as  a  function  of  distance  along  the  antenna 


is  then  computed  by 
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I  (z)  =  2k  a  H  4  ia,z  ) 


for  d  <z  <h 


=  2  n  a 


n  =0 


vnanH 


(2).  . 

i  (vna)cos 


m  =  0 


H  o  (vm  a  ) 


{Eq.  1-31} 

The  input  impedance  is  then  given  as  the  ratio  of  the  gap 
voltage  Vo  to  the  current  at  the  end  of  the  gap 


Z.  -  — Yp. 

I  (2  =d) 


{Eq.  1-32} 


B.  SOFTWARE  OVERVIEW 

The  computer  program  was  written  in  FORTRAN  77  and  was 
designed  to  run  on  an  MS-DOS  compatible  system  with  a  numeric 
coprocessor  and  460  kilobytes  of  available  memory  after  the 
disk  operating  system  in  loaded.  If  the  numeric  coprocessor 
is  not  available,  the  source  code  could  be  recompiled  to 
emulate  the  specific  coprocessor  calls.  The  user  would  notice 
a  dramatic  increase  in  run  time,  but  as  all  calculations  are 
done  in  single  precision,  the  accuracy  should  not  be  affected. 

The  program  code  is  divided  into  two  separate  programs. 
The  first  program,  called  MONO. FOR,  is  the  main  program  that 
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calls  a  number  of  subroutines  in  SUBS. FOR.  Variable  names 
used  in  the  program  code  closely  follow  those  used  in  the 
theoretical  development.  Each  variable  is  explicitly  declared 
in  each  subroutine  and  array  sizes  are  kept  to  a  minimum  to 
reduce  required  memory  at  run  time. 

The  first  section  of  code  in  MONO. FOR  is  used  to  input 
the  physical  parameters  using  either  an  unformatted  ASCII  data 
file  as  shown  in  Appendix  E  or  by  entering  the  values  from  the 
keyboard.  In  each  case,  measurements  may  be  entered  in  either 
wavelength  or  meters.  If  a  data  file  is  used,  the  first  line 
must  be  either  a  "w"  or  Hm" ,  and  may  be  either  upper  or  lower 
case,  to  distinguish  between  wavelengths  or  meters.  Unit 
number  zero  was  used  for  all  READ  statements,  which  allows  use 
of  the  program  in  a  batch  mode  where  all  inputs  and  screen 
directed  outputs  can  be  redirected  by  standard  DOS  methods. 
However,  the  use  of  compilers  other  than  Microsoft  Fortran 
(Versions  4.01  or  4.1)  may  produce  compiler  errors  when  using 
this  designation.  The  user  should  consult  the  specific 
compiler  manual  for  details  on  unit  designations. 

Before  prompting  the  user  to  identify  the  output  file 
name,  the  program  checks  for  two  conditions  that  will  produce 
inaccurate  results.  The  first  condition  is  when  the  upper 
plate  height  is  a  multiple  of  one  half  wavelength.  Physical¬ 
ly,  this  condition  will  allow  resonant  cavity  modes  between 
the  plates,  introducing  instabilities  in  the  numerical 
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solution.  The  mathematical  basis  for  this  can  be  easily 
shown.  When  the  upper  plate  height  is  a  multiple  of  half  the 
wavelength,  i /„  will  become  zero  for  certain  values  of  n.  As 
can  be  seen  in  Table  I,  the  value  of  a  Hankel  function  ap¬ 
proaches  - joo  as  the  argument  approaches  zero.  This  test 
occurs  in  the  main  program  and  supplies  earlier  error  detec¬ 
tion.  The  singularity  is  not  checked  for  in  the  subroutine 
that  computes  the  Hankel  functions  (HAN1)  found  in  SUBS. FOR. 

The  second  condition  that  is  tested  is  when  the  distance 
of  the  top  ground  plane  above  the  antenna,  designated  q,  is  an 
integer  multiple,  of  the  antenna  height.  This  is  similar  to 
the  previous  case  in  that  a  resonant  area  is  present  between 
the  antenna  and  the  upper  ground  plane.  However,  in  this 
case,  the  value  of  u. n,  which  appears  only  in  the  k  matrix, 
becomes  zero  for  multiples  of  the  ratio  h  to  q.  Figure  3 
illustrates  the  behavior  of  the  Bessel  function  for  a  range 
of  arguments  [Ref.  4],  The  zero  order  Bessel  function,  J0, 
is  finite  for  small  arguments;  however  the  first  order  Bessel 
function,  Ju  (found  in  the  denominator  of  Equation  1-30) 
rapidly  approaches  zero  as  the  argument  approaches  zero, 
causing  the  errors. 

In  an  effort  to  reduce  computation  time  at  the  expense  of 
memory  requirements,  several  arrays  are  dynamically  loaded  so 
that  the  values  may  be  easily  extracted  to  fill  the  system 
matrix  A* >k  and  the  driving  vector  Bb.  These  include  the 
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Figure  3  Bessel  Functions  of  the  First  Kind 


values  for  the  Fourier  moment  integrations,  the  Bessel 
function  values  and  the  Pm<k  matrix.  The  T  matrix  has  dimen¬ 
sions  of  x  Nx  and  consists  of  terms  from  Equation  1-23. 
The  I  matrix  has  dimensions  of  Mx  x  3  where  the  column 
positions  correspond  to  the  three  equations  in  Equation  1-18. 
The  Ja  and  Pmk  matrices  are  loaded  in  the  subroutine  ALOAD  and 
passed  back  to  the  main  program.  The  Ja  matrix  is  Mx  x  4  in 
size  where  the  column  positions  represent  the  following 
values: 
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1 


Position 

2 

3 

4 


Value 
Uj,*  JofUn  a) 

A  (Un  a) 

vn*  H1  (vn  O) 

( vj 2  *  H0  (Vtl  a) 


For  large  values  of  Nj,,  N2  and  Mx,  about  60%  of  the  total 
computation  time  is  involved  in  filling  the  last  temporary 
array,  the  P„,k  matrix.  This  array  is  dimensioned  to  Nx  x 
and  uses  values  from  all  three  arrays  for  each  location.  By 
filling  these  arrays  first,  the  computation  time  was  decreased 
by  an  average  of  85%,  however  the  memory  required  increased 
by  300%! 

Once  the  system  matrix  and  driving  vector  are  stored  in 
memory,  it  is  a  simple  task  to  solve  for  the  unknown  a^* s 
utilizing  the  back  substitution  method.  First  the  matrix  is 
upper  triangulated  by  performing  an  L-U  decomposition,  with 
pivoting,  on  the  An  t  matrix  using  the  FACTOR  subroutine.  Once 
this  is  accomplished,  the  SOLVE  subroutine  calculates  the 
eigenvalues  by  back  substitution. 

To  maintain  the  convention  that  an  outwardly  travelling 
wave  has  a  negative  imaginary  exponential  argument,  the  CSR 
function  was  included  in  SUBS. FOR.  since  not  all  compilers 
would  return  a  negative  imaginary  value  from  the  square  root 
algorithm  when  the  argument  is  a  negative  real  value,  the  CSR 
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function  insures  the  -j  component  of  the  square  root  of  a 
negative  argument  is  returned. 

The  subroutines  that  compute  the  Bessel  and  Hankel 
functions  for  complex  arguments  and  the  matrix  solution 
algorithms  were  supplied  by  Professor  Michael  A.  Morgan  of  the 
Naval  Postgraduate  School.  The  Bessel  and  Hankel  subroutines 
utilize  the  direct  power  series  method  for  small  values  (less 
than  or  equal  to  five  for  the  Bessel  subroutine  and  six  for 
the  Hankel  subroutine)  and  Hankel' s  asymptotic  formula  for 
larger  values. 

The  output  is  directed  to  a  user  specified  file  on  any 
drive  or  subdirectory  as  long  as  the  entire  string  is  less 
than  25  characters.  Longer  strings  will  truncate  the  extra 
characters  without  warning.  Printed  output  includes  all  input 
variables  including  the  number  of  coefficients  chosen  for 
Regions  I  and  II  and  for  the  gap  voltage  expansion.  Length 
measurements  are  reported  in  both  meters  and  fractions  of  a 
wavelength.  The  system  eigenvalues  are  then  printed  followed 
by  the  position  number,  distance  from  the  end  of  the  gap  and 
magnitude  and  phase  of  the  computed  currents,  up  to  the  upper 
plate  height.  Finally  the  input  resistance  and  reactance  is 
computed  and  printed  to  the  file.  Sample  input  and  output 
files  are  listed  in  Appendix  E. 
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III.  VALIDATION 


In  order  to  validate  any  computational  algorithm,  one 
must  establish  the  accuracy  obtained  by  comparing  the  computed 
results  against  either  experimental  data  or  another  widely 
accepted  software  package  that  has  been  proven  to  be  accurate. 
For  this  case,  the  logical  choice  was  to  use  the  Numerical 
Electromagnetics  Code  to  calculate  the  input  impedance  and 
current  distribution  for  a  range  of  different  physical 
structures  and  compare  these  results  against  those  obtained 
by  MONO.  From  this  point,  the  program  sensitivity  will  be 
examined  for  various  ground  plate  spacings,  antenna  lengths 
and  radii.  Chapter  4  will  investigate  the  effect  that 
reducing  the  number  of  modal  expansion  terms  has  on  this 
accuracy . 

A.  THE  NUMERICAL  ELECTROMAGNETICS  CODE 

The  Numerical  Electromagnetics  Code  (NEC)  is  widely 
accepted  as  an  accurate  method  of  computing  the  antenna 
currents,  input  impedance  and  radiation  patterns  of  antennas 
that  can  be  modeled  using  wire  or  patch  structures.  NEC 
calculates  these  values  by  solving  integral  equations  for  the 
currents  imposed  on  the  structure  by  a  voltage  source  or  an 
incident  electromagnetic  wave.  The  structure  is  best  modeled 
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by  straight  wire  segments  that  should  be  less  than  ten  percent 
of  the  operating  wavelength  and  conforms  to  the  physical 
geometry  using  a  best  fit  approach.  [Ref.  6,  p.  3] 

The  user  supplies  NEC  with  a  formatted  data  file  that 
specifies  the  location  of  the  structure  in  the  rectangular 
coordinate  system,  the  number  of  segments  used  to  compute  the 
results,  the  driving  source  location  and  output  parameters  and 
the  existence  and  characteristics  of  a  ground  plane  in  the  XY 
plane.  Because  of  the  infinite  number  of  images  produced  by 
a  pair  of  perfect  ground  planes,  the  user  must  also  construct 
enough  images  with  the  proper  driving  voltages  that  will  yield 
consistent  results.  NEC  then  calculates  the  segment  length, 
and  uses  this  value  as  the  gap  distance  between  the  ground 
plane  and  the  base  of  the  antenna.  The  currents  are  calcu¬ 
lated  at  the  midpoint  of  each  segment  vice  at  the  ends  as  is 
done  in  MONO. 

NEC  uses  two  approximations  in  the  calculation  of  the 
electric  field,  the  thin-wire  kernel  and  the  extended  thin- 
wire  kernel.  The  thin-wire  kernel  can  be  used  when  the 
segment  length  to  radius  ratio  is  greater  than  8,  while  the 
extended  wire  kernel  is  accurate  when  this  ratio  is  at  least 
2.  Since  NEC  sets  the  gap  distance  equal  to  the  segment 
length  and  accuracy  is  improved  if  the  segment  lengths 
(especially  across  a  boundary)  are  of  equal  size,  setting  the 
gap  distance  equal  to  the  antenna  diameter  forces  the  segment 
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length-to-radius  ratio  to  be  equal  to  2.  Additionally,  this 
gap  distance  is  consistent  with  the  requirement  for  a  constant 
electric  field  potential  inside  the  gap  region.  Therefore, 
the  extended  wire  kernel,  which  uses  a  uniform  surface  current 
along  the  segment  length  and  assumes  no  variation  of  the 
currents  along  the  <p  direction,  will  be  used  for  all  valida¬ 
tion  runs.  These  approximations  also  support  the  basic 
assumption  of  the  TM  mode  dominating  the  currents  on  the 
monopole.  For  the  conditions  listed  above,  NEC's  accuracy  is 
within  1%  of  experimental  results.  [Ref.  6] 

B.  COMPARISONS  MADE 

The  first  goal  will  be  to  determine  the  minimum  number  of 
reflections  that  must  be  assembled  in  the  NEC  data  set  to 
obtain  consistent  results  for  a  quarter  wave  antenna  struc¬ 
ture.  From  this  data,  an  appropriate  number  of  reflections 
will  be  used  that  represents  consistent  results  while  reducing 
the  calculation  time  for  all  future  runs.  This  will  be 
followed  by  demonstrating  consistent  results  for  various 
ground  plane  spacings.  The  upper  ground  plane  was  introduced 
to  establish  a  closed  region.  The  upper  plate  should,  if  far 
enough  away  from  the  antenna,  have  little  effect  on  the  input 
impedance  and  antenna  currents  for  distances  that  are  not 
exact  multiples  of  half  the  wavelength.  Again  a  value  of  / 
will  be  selected  that  represents  reliable  results  but  reduces 


27 


the  calculation  time  required  for  MONO.  The  antenna  radius 
will  then  be  adjusted  over  the  range  of  0.005A  to  0.02 a 
followed  by  variations  in  the  antenna  height  for  a  selected 
radius . 

For  all  MONO  calculations,  the  number  of  modes  used  in 
Regions  I  and  II  (N1#  N2)  will  always  be  set  to  60,  while  the 
number  of  modal  expansion  terms  used  for  the  gap  field  (Mx) 
will  be  set  to  two  times  the  upper  plate  spacing  to  gap 
distance  ratio.  Chapter  4  will  further  investigate  the 
accuracy  of  the  results  as  a  function  of  the  number  of  modes 
used.  A  preliminary  check  for  consistent  results  can  be 
obtained  by  evaluating  the  differences  in  the  input  impedance 
results.  But  these  values  may  differ  slightly  since  NEC 
calculates  the  current  at  the  midpoint  of  each  segment  instead 
of  at  the  ends  of  the  segment  as  is  done  in  MONO.  Addition¬ 
ally,  since  the  first  segment  used  in  NEC  is  the  gap  distance, 
the  input  impedance  is  actually  calculated  for  a  point  not  on 
the  monopole  but  in  the  center  of  the  feed  region.  In 
comparison,  MONO  calculates  the  input  impedance  at  the  base 
of  the  antenna.  The  effect  this  has  on  the  comparisons  will 
be  discussed  later  in  this  chapter.  Additionally,  these 
differences  may  be  emphasized  since  the  magnitude  of  the 
current  at  the  end  of  the  gap  region  can  be  quite  small  so 
that  even  small  differences  in  the  calculated  values  will 
yield  relatively  high  differences.  Therefore,  to  better 
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assess  the  accuracy  of  the  results,  the  calculated  current 
distribution  along  the  antenna  will  be  compared  graphically 
with  that  obtained  from  NEC. 

C.  RESULTS 

Since  the  ground  planes  used  to  divide  the  area  into 
regions  create  an  infinite  number  of  reflections,  one  must 
first  examine  the  minimum  number  of  reflections  required  to 
obtain  convergence  to  the  half  space  case,  having  no  upper 
ground  plane.  For  this  set  of  runs,  the  antenna  is  0.24*  high 
with  radius  of  0.01A  and  the  upper  ground  plate  was  placed  1.4 
wavelengths  from  the  lower  plate.  Appendix  F  has  a  complete 
listing  of  some  of  the  NEC  data  sets  used  for  the  validations. 
It  should  be  noted  that  the  distance  between  the  lower  plate 
and  the  center  of  the  reflected  dipole  is  equal  to  twice  the 
plate  spacing.  This  is  due  to  the  complete  imaging,  including 
the  lower  ground  plane,  of  the  physical  structure  from  the 
upper  ground  plane.  One  can  see  from  Figure  4  that  adding  a 
second  reflection  does  not  change  the  current  distribution  by 
any  appreciable  amount,  therefore  one  reflection  will  be  used 
for  all  further  validations. 

Now  that  a  baseline  data  set  has  been  decided  upon,  the 
effects  of  various  plate  spacing,  antenna  heights  and  radius 
will  be  examined.  The  first  item  to  look  at  will  be  the 
effect  the  plate  spacing  has  on  the  current  distribution  and 
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Antenna  Current  Distribution  (NEC) 
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Figure  4  Antenna  Current  vs  Number  of  Reflections 
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the  calculated  input  impedance.  Table  II  shows  the  numerical 
results  for  an  antenna  of  height  0.24 A  and  radius  of  0.01A 
placed  between  ground  planes  that  are  spaced  3.1,  1.4  and  0.81 

Table  II  RESULTS  FOR  VARIOUS  UPPER  PLATE  SPACINGS 


Plate 

Spacing 

ES 

mam 

MOl 

Rin 

JO 

Xin 

m 

a 

3.10  A 
1.40  A 
0.81  A 

44.63 

44.90 

43.51 

13.35 

13.97 

14.84 

42.65 

43.14 

44.67 

13.14 

13.95 

13.84 

4.64 

4.08 

2.66 

m 

wavelengths.  Figures  5  through  7  show  the  current  distribu¬ 
tion  for  each  of  the  above  cases.  While  none  of  the  results 
exceed  acceptable  limits,  the  case  of  /=1.4A  is  closer  to  the 
results  computed  by  NEC  and  will  be  used  as  the  baseline  for 
further  comparisons.  It  should  be  noted  here  that  while  the 
example  values  selected  for  the  plate  spacing  have  little 
effect  on  the  result,  the  plate  spacing  can  not  be  an  integer 
multiple  of  a  half  wavelength.  Likewise,  the  distance  between 
the  top  of  the  antenna  structure  and  the  upper  plate  can  not 
be  a  multiple  of  a  half  wavelength.  Physically,  these 
configurations  would  present  the  ability  to  support  resonant 
modes  that  do  not  require  a  driving  source.  In  reality,  this 
condition  will  not  exist  as  there  is  always  a  source  of 
resistance  which  causes  all  modes  to  decay  in  time.  However, 
the  first  condition  will  exist  numerically  where  H0<2)  ap¬ 
proaches  negative  infinity  when  vB  (Equation  1-11)  equals 
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NEC  vs  MONO 
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zero.  By  replacing  the  wave  number  k0  with  2n/x,  then  is 
zero  when  the  plate  spacing  /  is  an  integer  multiple  of  a  half 
wavelength.  A  similar  condition  will  also  occur  when  the 
distance  between  the  top  of  the  structure  and  the  upper  plate 
is  a  multiple  of  a  half  wavelength.  In  this  case,  Jx  found  in 
the  denominator  of  PB  k  (Equation  1-30)  is  zero  when  is  zero. 
The  program  code  will  alert  the  user  to  both  conditions  and 
halt  execution. 

The  next  variable  to  consider  is  the  antenna  height.  For 
extremely  short  antenna  structures,  the  current  decays  almost 
linearly  along  the  wire  from  its  initial  value  to  zero.  When 
the  antenna  is  about  one  quarter  wavelength  long,  the  current 
distribution  takes  on  the  familiar  cosine  shape.  As  the 
antenna  length  approaches  one  half  wavelength,  the  distribu¬ 
tion  approaches  a  sine  wave  shape  and  the  input  impedance 
approaches  infinity  as  the  current  approaches  zero. 

Table  III  lists  the  results  when  the  structure  is  placed 
between  the  two  ground  planes  spaced  1.4 A  apart  for  various 
antenna  heights  and  radii.  In  the  case  where  the  antenna 
height  approaches  a  half  wavelength,  the  percent  difference 
in  the  input  impedance  values  appears  to  be  large.  This  can 
be  attributed  to  the  differences  in  the  location  on  the 
structure  where  the  input  impedance  is  actually  calculated. 
This  disparity  combined  with  the  low  magnitudes  in  general  can 
lead  to  alarming  discrepancies  with  the  NEC  results. 
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Table  III  NEC,  MONO  RESULT  COMPARISONS  FOR  VARIOUS  h  AND  a 


mm 

Mono 

■E 

9E 

mEBm 

test 

Rin 

Xin 

mm 

Mm 

.12 

.005 

6.6 

-137.9 

7.4 

-148.7 

10.7 

.24 

.005 

42.0 

12.4 

40.6 

10.8 

3.3 

.48 

.005 

254.8 

-285.0 

339.9 

-281.7 

25.0 

KQ 

.12 

.01 

1 

-104.6 

9.3 

-118.2 

17.8 

11.5 

.24 

.01 

14.0 

43.1 

14.0 

4.0 

0.1 

.48 

.01 

-203.8 

264.8 

-199.9 

32.9 

2.0 

.12 

.02 

9.8 

-76.3 

13.4 

-92.2 

.24 

.02 

47.9 

13.4 

46.1 

15.9 

.48 

.02 

129.6 

-139.8 

220.1 

-120.9 

41.1 

1  :  1 

Therefore,  one  must  look  at  the  current  distribution  to 
accurately  determine  the  differences  in  the  computed  results. 
Figures  8  through  11  show  that  the  differences  in  the 
magnitude  of  the  current  distribution  are  within  two  percent. 
Additionally,  one  could  also  compare  the  phase  difference 
between  the  two  methods  to  further  establish  an  acceptable 
method  of  judging  the  accuracy.  Figure  12  shows  that  the 
difference  between  the  phase  of  the  current  as  calculated  by 
MONO  and  NEC  is  well  within  two  percent  along  the  entire 
length  of  the  antenna. 

Another  point  that  can  be  observed  is  the  effect  that  a 
change  in  the  antenna  radius  has  on  the  results.  Once  again 
a  review  of  the  current  distribution  is  required  to  accurately 
assess  the  effect  any  errors  would  have  on  the  far  field 
pattern;  however  one  detail  is  evident.  As  the  antenna 
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becomes  thinner  the  difference  between  the  two  results  is 
reduced.  This  effect  can  be  attributed  to  the  increased 
resolution  in  the  NEC  results  as  the  radii  and,  therefore,  the 
gap  distance  are  reduced.  Additionally,  the  approximation  to 
a  TM  current  distribution  along  the  surface  of  the  structure 
becomes  more  accurate  with  a  thinner  structure. 

In  review,  it  has  been  shown  that  due  to  the  differences 
in  the  methods,  the  best  method  to  assess  the  accuracy  of  MONO 
using  NEC  is  obtained  by  comparing  both  the  current  magnitude 
and  phase  distribution  along  the  antenna  surface  vice  the 
input  impedance  values.  Using  this  criterion,  MONO  is 
consistently  within  two  percent  of  the  results  obtained  with 
NEC  for  a  variety  of  structures.  Chapter  4  will  investigate 
the  difference  caused  by  reducing  the  number  of  expansion 
terms  and  attempt  to  identify  a  minimum  number  of  terms 
required  to  obtain  results  that  are  consistent  with  those 
obtained  here. 
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Figure  8  NEC  vs  MONO  for  Quarter  Wave  Monopole 
(Radius  =  0.005A) 
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IV.  CONVERGENCE  CONSIDERATIONS 


In  any  realizable  computing  machine  where  there  is  a 
finite  quantity  of  storage,  the  first  question  that  must  be 
answered  is  the  number  of  terms  required  to  obtain  accurate 
results.  If  one  uses  too  few  terms,  the  results  are  er¬ 
roneous.  However,  using  too  many  terms  can  often  lead  to 
inaccuracies  as  well.  If  the  additional  terms  are  small 
enough  to  be  approximately  on  the  same  order  as  the  precision 
of  the  computer,  then  the  additional  terms  may  actually 
increase  the  error.  To  compound  the  problem,  the  additional 
terms  require  additional  memory  and  increase  the  time  required 
to  compute  the  results.  In  investigating  this  question,  we 
will  first  examine  the  number  of  terms  required  for  Regions 
I  and  II  (Nx  and  N2)  followed  by  an  analysis  of  the  number  of 
terms  required  to  represent  the  electric  field  in  the  gap 
region  (Mx) . 

A.  MODAL  TRUNCATIONS  FOR  REGIONS  I  AND  II 

In  the  previous  chapter,  the  number  of  modes  chosen  for 
the  two  regions  was  the  maximum  allowed  in  the  program,  namely 
60.  But  is  that  enough?  One  could  say  that  since  the  results 
compare  favorably  with  NEC,  the  solution  has  converged.  A 
more  accurate  approach  is  to  increase  the  number  of  modes  and 
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compare  the  results.  However,  due  to  the  memory  limit  of  640 
kilobytes  of  addressable  memory  imposed  by  the  compiler,  the 
source  code  was  recompiled  on  a  80386-specific  FORTRAN 
compiler  that  allows  for  larger  arrays  and  is  only  limited  by 
the  available  memory  of  the  computer  being  used.  Figure  13 
compares  the  current  distribution  of  a  quarter  wave  antenna 
that  has  the  upper  plate  height  set  at  a  large  value  (3.1A) 
for  60  modes  and  100  modes.  This  additional  information 
proves  that  the  system  has  converged  to  an  accurate  solution 
with  60  modes  and  any  additional  modes  are  not  required.  The 
effect  reducing  this  number  has  on  the  accuracy  will  be 
discussed  later  in  this  chapter. 

B.  MODAL  NEEDS  FOR  THE  GAP  REGION 

As  mentioned  earlier,  the  original  approach  was  to  expand 
the  electric  field  along  the  surface  of  the  antenna  in  one 
expansion  between  the  ground  planes  as  shown  below 


1S 
d  ' 


for  0  <z  <d 


0, 


for  d£z  <h 


l  X  uncn  '<,(“„<*)  cos  nK^a  —  ,  for  hZzZl 
v  H 


0 


{Eq.  4-1} 
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However,  this  approach  failed  to  yield  consistent  results  and 
a  review  of  the  source  code  showed  no  errors  in  the  logical 
design.  Upon  an  examination  of  the  electric  field  in  Region 
I,  it  was  discovered  that  a  very  high  value  for  Ni  would  be 
required  in  order  to  obtain  convergence  to  the  known  field  in 
the  narrow  feed  gap.  It  was  believed  that  significant  errors 
in  the  electric  field  would  cause  inaccuracies  in  the  result. 
The  errors  in  the  results  were  thought  to  be  caused  by  the 
lack  of  terms  used  to  describe  the  electric  field.  The 
electric  field  in  the  gap  region  can  be  thought  of  as  a  step 
function  of  short  duration.  In  the  frequency  domain,  func¬ 
tions  of  this  type  occupy  a  large  spectral  bandwidth,  requir¬ 
ing  many  terms  for  the  description.  The  large  expected  value 
of  Ni  would  necessitate  the  inversion  of  a  very  large  matrix, 
thus  making  this  approach  inefficient.  To  avoid  this  problem, 
a  separate  expansion  for  the  constant  electric  field  was  added 
to  Equation  4-1  to  form  Equation  1-14. 

Once  again  the  number  of  terms  required  was  the  main 
question.  However,  this  field  is  known  prior  to  the  computa¬ 
tion  and  the  number  of  terms  can  be  estimated  using  standard 
Fourier  analysis.  The  Fourier  transform  of  a  step  function 
is  known  to  be  a  sine  function  with  zero  crossings  at  multi¬ 
ples  of  1/d,  where  d  is  the  duration  of  the  pulse.  Each  lobe 
consists  of  impulse  functions  at  intervals  of  the  fundamental 
frequency,  /„,  which  equals  1/T0,  where  T0  is  the  period  of  the 


46 


pulse  train.  In  this  case,  the  upper  ground  plane  causes  a 
reflection  of  the  original  structure,  but  with  opposite 
polarity,  at  a  distance  of  21  from  the  lower  ground  plane. 
Therefore  the  number  of  impulses  in  each  lobe  of  the  sine 
function  equals  2 l/d.  This  can  be  investigated  graphically  as 
well.  Figures  14  to  19  (pages  48  to  50)  illustrate  the  field 
computed  for  a  structure  with  a  gap  distance  ( d )  of  0.06 
meters  and  upper  plate  height  of  9.54  meters  and  fed  by  a 
constant  one  volt  source  for  various  numbers  of  terms  in  the 
series.  As  can  be  expected,  when  the  number  of  terms  are  not 
an  integral  multiple  of  l/d,  excess  oscillations  occur.  It  is 
also  important  to  note  that  as  the  number  of  terms  increase, 
the  calculated  curve  approaches  the  ideal  value  at  the 
discontinuities,  while  the  pointwise  error  along  the  rest  of 
the  curve  increases.  However,  the  error  in  the  least  square 
sense  decreases.  This  is  known  as  the  Gibbs  phenomenon.  Near 
a  discontinuity  for  a  large  number  of  terms,  the  overshoot  is 
approximately  8.95  percent  of  the  desired  value,  or  -18.17 
volts  (Ref.  5  p.  557].  This  value  is  reached  when  the  number 
of  terms  equals  3 l/d  as  seen  in  Figure  19,  therefore  addition¬ 
al  terms  will  not  reduce  the  error  at  the  discontinuities. 
From  these  graphs,  one  can  see  that  the  number  of  terms  used 
to  represent  the  electric  field  in  the  gap  region  should  be 
an  integer  multiple,  less  than  or  equal  to  three,  of  the  ratio 
of  upper  plate  height  to  gap  distance.  An  upper  plate  height 
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E  Field  Convergence 

For  100  modes  (0.631/d) 


Figure  14  Gap  Field  Representation  for  100  Terms 


E  Field  Convergence 
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Figure  15  Gap  Field  Representation  for  159  Terms 
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Figure  17  Gap  Field  Representation  for  318  Terms 
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Figure  18  Gap  Field  Representation  for  400  Terms 
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Figure  19  Gap  Field  Representation  for  477  Terms 


to  gap  distance  ratio  of  250  was  considered  to  be  the  maximum 
desired.  Therefore  to  insure  the  code  would  fit  easily  into 
the  640  kilobyte  limit,  the  number  of  terms  used  for  this 
expansion  was  set  at  21/ d.  This  was  the  value  used  for  all 
the  validation  calculations  completed  in  Chapter  3. 

C.  EFFECT  OF  REDUCING  N1  AND  N2 

If  the  number  of  terms  can  be  reduced  and  still  provide 
consistent  results,  one  can  then  decrease  the  size  of  the 
internal  arrays  and  reduce  the  computation  time.  For  example, 
by  reducing  Nx  by  50%,  the  computation  time  is  reduced  by 
almost  48%,  and  when  both  Nj  and  N2  are  reduced  by  50%,  the 
computation  time  is  reduced  by  71.5%!  However,  the  current 
distribution  along  the  antenna  is  unknown,  making  it  difficult 
to  calculate  the  number  of  coefficients  required  for  consist¬ 
ent  results  based  on  the  input  parameters.  The  results 
obtained  by  independently  varying  Nx  and  N2  are  shown  in  Table 
IV  and  are  compared  to  the  results  obtained  in  Chapter  3  for 
an  antenna  with  the  following  physical  parameters: 


Antenna  Height  (h) 

0 . 24  A 

Antenna  Radius  (a) 

0.01A 

Gap  Distance  ( d ) 

0 . 02  A 

Plate  Spacing  (/) 

2 . 17  A 

Number  of  terms  for 
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Gap  Region  (Mj) 
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Tabl*  IV  COMPARISON  OF  RESULTS  AS  Nx  AND  N2  VARY 


N1 

N2 

Rin 

Xin 

■B 18 

60 

60 

43.08 

13.98 

N/A 

N/A 

50 

50 

42.70 

13.50 

0.88 

3.43 

40 

40 

42.80 

13.10 

0.65 

6.29 

30 

30 

41.60 

12.40 

3.44 

11.30 

20 

20 

41.20 

9.80 

4.36 

29.90 

50 

60 

42.20 

12.50 

2.04 

10.59 

40 

50 

41.80 

11.30 

2.97 

19.17 

30 

40 

39.90 

8.60 

7.38 

38.48 

20 

30 

35.60 

-0.80 

17.36 

105.72 

60 

50 

44.60 

16.70 

3.53 

19.46 

50 

40 

44.90 

17.47 

4.22 

24.96 

40 

30 

45.70 

18.99 

6.08 

35.84 

30 

20 

47.40 

22.60 

10.03 

61.66 

As  expected,  the  difference  in  the  input  impedance  increases 
as  the  number  of  terms  is  reduced.  However,  the  interesting 
point  is  the  effect  that  the  apparent  relationship  between  Nx 
and  Nz  has  on  the  result.  If  Nj  and  N2  are  equal,  acceptable 
results  are  obtained  for  as  low  as  30  terms,  while  one  must 
use  at  least  50  terms  if  the  values  are  not  equal,  to  get 
comparable  results.  A  mathematical  explanation  for  this 
relationship  is  not  yet  understood  and  will  be  investigated 
in  future  thesis  efforts.  Figure  20  represents  the  current 
distribution  as  N2  and  N2  vary  independently  compared  with  the 
result  when  Ni  and  N2  equal  60. 


D.  EFFECT  ON  ACCURACY  AS  M,  VARIES 

As  seen  in  the  previous  section,  one  can  reduce  the 
memory  requirements  and  computation  time  by  reducing  the 
number  of  coefficients  used  in  Regions  I  and  II.  Using 
Fourier  analysis,  the  user  can  calculate  a  value  for  the 
number  of  terms  required  to  accurately  represent  the  electric 
field  along  the  gap  (Mx) .  Since  this  value  is  much  larger 
than  Nx  and  N2,  one  could  dramatically  reduce  the  required 
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storage  and  speed  up  the  computations  by  reducing  Mx  to  the 
same  order  as  Nt  and  N2. 


To  demonstrate  the  sensitivity  of  the  input  impedance  to 
the  number  of  terms  used  for  the  electric  field  along  the  gap, 
one  can  compare  the  results  against  those  found  when  Nx  and  Nz 
are  60  terms  and  Hx  is  computed  as  21/ d.  Table  V  lists  the 
results  for  same  antenna  structure  described  on  page  51.  As 
expected,  the  errors  in  the  calculated  input  impedance  begin 
to  increase  as  the  number  of  terms  are  reduced.  However,  this 
difference  is  not  linear  for  all  values  of  Nx  and  Nz.  In  fact, 
Figure  21  suggests  an  optimal  number  of  Mx  terms,  for  given 
Nz  and  Nz,  exists  which  results  in  smaller  differences  between 
the  two  cases.  Since  is  a  function  of  the  gap  distance, 


Table  V  MONO  SENSITIVITY  TO  (Mx)  FOR  GIVEN  Nx  AND  Nz 


N1 

N2 

Ml 

R 

X 

mm 

mm 

60 

60 

ESI 

43.08 

13.98 

BBS 

N/A 

60 

60 

9 B3 

43.31 

13.87 

i 

0.77 

60 

60 

El 

43.02 

13.85 

0.92 

60 

60 

40 

43.39 

13.95 

0.72 

0.20 

60 

60 

30 

42.44 

14.14 

1.49 

1.16 

50 

50 

60 

42.81 

13.44 

0.63 

3.85 

50 

50 

50 

42.59 

13.42 

1.14 

3.99 

50 

50 

40 

43.12 

13.53 

0.09 

3.21 

50 

50 

30 

42.13 

13.72 

2.21 

1.85 

30 

30 

40 

41.72 

12.39 

3.16 

11.36 

30 

30 

30 

40.70 

12.64 

5.52 

9.57 

30 

30 

20 

42.02 

13.26 

2.46 

5.14 

30 

30 

15 

39.41 

13.61 

8.52 

2.63 

30 

30 

10 

38.29 

14.73 

11.12 

5.38 
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the  input  impedances  were  recalculated  when  the  gap  distance 
is  reduced  to  0.01A  for  selected  values  of  N*  and  N2  as  shown 
in  Figure  22.  While  this  clearly  demonstrates  that  an  optimal 
number  of  terms  may  exist  that  produces  results  which  approach 
those  for  large  values  of  Mu  Ni  and  N2,  a  clear  relationship 
between  the  physical  structure  and  these  values  are  not  yet 
fully  understood. 

But  how  can  the  number  of  terms  be  reduced  by  almost  91 
percent,  and  yet  the  difference  is  no  greater  than  5.13 
percent?  First,  the  basis  for  selecting  the  number  of  terms 
used  was  based  on  obtaining  accurate  representation  for  the 
electric  field  along  the  antenna.  However,  since  the  far 
field  patterns  are  derived  directly  from  the  currents,  it  is 
the  current  that  must  remain  consistent  and  not  necessarily 
the  electric  field  along  the  surface  of  the  antenna.  From 
Equation  1-31,  we  see  that  the  current  is  derived  from  the  H^, 
or  the  first  derivative  of  the  potential  field  with  respect 
to  p.  However,  the  electric  field  along  the  z  direction  is  a 
function  of  the  second  derivative  of  the  potential  with 
respect  to  z.  Since  the  derivative  process  increases  the 
noise,  or  errors,  it  is  evident  that  small  differences  in  the 
potential  field  may  yield  wide  variations  in  the  electric 
field  while  the  magnetic  field,  and  therefore  the  current 
distribution,  may  remain  within  accepted  accuracy  limits.  To 
prove  this  point,  Equation  4-2  shows  an  expression  that 
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Average  Error  (percent)  J  |  Average  Error  (percent) 


Average  Error  vs  Ml 

For  given  number  of  terms  N1  <=  N2) 


N 1  =60 
N 1  =40 
N1=30 


10  20  30  40  50  60 

Number  of  Modes  Used  for  Gap  Region  (Ml) 

Figure  22  Average  Differences  for  Reduced  Values  of  Mi  for 
Given  N2  and  N2.  (Gap  =  0.01 A) 


approximates  the  electric  field  along  the  surface  of  the 
antenna  by  taking  the  second  derivative  of  Equation  1-31  with 
respect  to  z. 


A/, 


V'o 

+  JQ)£  OjZ 


m  =  0 


m 


V. 


{Eq.  4-2} 

This  equation  can  now  be  used  to  visualize  the  approximate 
change  in  the  electric  field  along  the  antenna  for  a  given 
number  of  modes  MA.  As  shown  in  Figure  23,  the  change  in  the 
magnitude  of  the  current  distribution  when  Mx  is  changed  from 
217  (which  equals  2l/d)  to  30  is  negligible.  However,  Figure 
24  shows  the  curves  obtained  from  Equation  4-2  for  the  same 
parameters  as  used  in  Figure  23.  As  one  can  see,  the  field 
is  accurately  represented  near  the  end  of  the  wire,  however 
near  the  gap  region,  there  are  an  insufficient  number  of 
coefficients  to  accurately  compute  the  fields  in  this  region. 
This  refutes  the  original  assumption  that  the  electric  fields 
must  be  accurately  represented  to  guarantee  an  accurate 
current  solution.  Additionally,  since  can  be  on  the  same 
order  as  the  values  used  for  Nx  and  N2,  then  the  coefficients 
used  to  describe  the  gap  field  and  the  field  along  the 
structure  can  be  combined  into  one  coefficient.  The  next 
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chapter  will  explore  the  theory  and  the  differences  this  new 
approach  has  over  calculating  a  separate  set  of  expansion 
coefficients  for  the  gap  region. 
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Figure  24  Approximated  Electric  Field  Distribution 


(Calculated  using  Equation  4-2) 
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V.  A  SIMPLER  APPROACH 


During  the  early  stages  of  the  software  development,  it 
was  decided  that  a  separate  expansion  for  the  constant 
electric  field  in  the  gap  region  was  required  to  obtain 
consistent  results  for  the  input  impedance  and  current 
distribution  along  the  structure.  This  decision  was  based  on 
the  belief  that  large  errors  in  the  total  electric  field  along 
the  structure  would  result  in  large  errors  in  the  current 
distribution.  However  since  the  current  distribution  is 
derived  directly  from  the  magnetic  field,  and  is  relatively 
insensitive  to  the  electric  field,  the  number  of  terms  used 
for  the  gap  region  expansion  was  reduced  to  the  same  order  as 
those  used  for  Nx  and  N2  without  introducing  appreciable 
errors.  Noting  this,  the  original  method  of  describing  the 
system  was  reevaluated  such  that  one  could  reduce  the  required 
memory  and  speed  up  the  calculations  even  further  by  combining 
the  a,,  and  bn  coefficients  used  in  Equation  1-10  (page  11)  into 
one  set  of  terms  and  by  deleting  the  separate  expansion  along 
the  gap  region  entirely.  This  chapter  will  work  backwards 
from  the  technique  presented  in  Chapter  2  and  will  note  the 
changes  in  the  system  equations  resulting  from  this  approach. 
The  resulting  code,  which  will  be  referred  to  as  M0N08,  will 
then  be  compared  to  the  results  obtained  for  the  cases  with 
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the  largest  differences  between  NEC  and  the  code  using  the 
separate  expansion  (called  M0N07) . 


A.  MODIFYING  THE  TECHNIQUE 

As  was  discussed  in  Chapter  2 ,  the  previous  approach  used 
two  expansions  in  Region  I  where  p  >  a.  The  first  expansion 
sets  the  tangential  electric  field  to  zero  from  the  lower 
ground  plane  to  the  top  of  the  antenna.  A  second  expansion 
sets  the  tangential  electric  field  to  -  V0/d  between  the  lower 
ground  plane  and  the  feed  point  of  the  antenna  and  sets  Ez  to 
zero  from  this  point  to  the  upper  ground  plane.  Since  the 
Fourier  transform  of  a  sum  is  the  sum  of  the  transforms  [Ref. 
7],  we  can  combine  the  expansions  for  the  potential  field  in 
Region  I  as  shown  below 


N, 


¥  (P  )  =  X  P  > C0S  ("7 


n  =  0 


{Eq.  5-1} 


where  a„  equals  the  sum  of  the  a„  and  bm  coefficients  used  in 
Chapter  2.  Again  we  match  the  magnetic  fields  across  the 
boundary  between  Region  I  and  II  such  that 


"x 


/M*>Z)=SX  VnanH\1)(y,P) 


COS 
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n  =o 
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{Eq.  5-2} 
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By  utilizing  the  orthogonality  property  of  the  Fourier  moment 
integrations  used  in  Equation  1-15,  the  system  equation  can 
be  expressed  ns 


N  i 

'4LivkakHf)(vka) 
1  =  0 


(3) 

m 


N  ■ 


In  =0 


0  fan  &  )  T  m„  T  kn 


-  vm2amH^\vma) 


{Eq.  5-3} 

for  m  =  0  to  Nx.  This  can  be  expressed  in  matrix  form  as 


Ny 

Am,fcai  =  for  m  =0,A/j 

1  =  0 

{Eq.  5-4} 

where  A„,_k  is  given  in  Equation  1-29.  As  the  fcm  coefficients 
are  now  included  in  the  an  terms,  the  driving  vector,  Bro,  is 
reduced  to 


B  =9/  (3)  — 

{Eq.  5-5) 

Once  the  unknown  an  coefficients  are  found,  the  current  along 
the  antenna,  as  a  function  of  z,  is  now  given  by 


N  , 

/  (2)  =2m  ^  vn  an H  {2)( vna)  cos  |-j-| 
n  =0 

{Eq.  5-6} 

B.  SOFTWARE  MODIFICATIONS 

As  seen  in  Equations  5-3  to  5-6  above,  removing  the 
separate  expansion  for  the  electric  field  along  the  gap  region 


62 


requires  changing  the  software  code  for  only  the  driving 
vector  Bm  and  the  expression  for  the  current  as  a  function  of 
position  along  the  antenna  structure.  The  first  advantage 
realized  was  the  reduction  in  the  amount  of  memory  storage 
required  to  run  the  program.  Since  it  was  shown  in  the 
previous  chapter  that  the  system  had  converged  when  the  number 
of  modes  in  Regions  I  and  II  were  set  at  60,  independent  of 
the  number  of  modes  in  the  gap  region,  it  was  not  necessary 
to  dimension  any  array  greater  than  60.  However,  we  can  now 
reduce  the  size  of  the  I,  T,  Ja,  and  P->k  matrices  by  almost  an 
order  of  magnitude,  such  that  the  amount  of  memory  required 
to  run  the  program  has  decreased  by  75  percent  to  less  than 
104  kilobytes.  Additionally,  the  execution  time  of  the 
program  has  decreased  by  an  average  of  35  percent  for  the  runs 
presented  here. 

But  does  this  method  yield  results  that  are  comparable  to 
the  results  when  the  separate  expansion  is  used?  We  will 
examine  the  differences  between  M0N07  and  M0N08  for  various 
plate  spacings,  antenna  heights,  as  well  as  the  effect  that 
independently  varying  N:  and  N2  has  on  the  results  obtained 
earlier. 

C.  RESULTS 

As  a  first  test  of  the  accuracy  for  M0N08,  we  will 
examine  the  differences  in  the  calculated  input  impedance 
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between  M0N07  and  M0N08  as  the  ground  plane  spacing  is  varied. 
The  three  structures  used  in  Chapter  3  for  Table  II  were  rerun 
where  Ni  and  N2  were  set  to  60.  Table  VI  indicates  that  the 


Table  VI  COMPARISON  OF  M0N07 ,  M0N08  FOR  VARIOUS  / 


Plate 

M0N07 

MONO  8 

Spacing 

Rin 

Xin 

Rin 

Xin 

3.10A 

42.65 

13.14 

43.09 

13.12 

1.03 

0.15 

1.40A 

43.14 

13.95 

43.15 

13.96 

0.02 

0.07 

0.81A 

44.67 

13.84 

44.68 

13.82 

0.02 

0.14 

results  obtained  from  the  new  method  differ  only  slightly  from 
the  results  obtained  when  the  gap  expansion  is  used. 

The  next  test  of  the  accuracy  considered  different 
antenna  heights.  In  this  case,  the  comparison  is  made  between 
the  structures  that  resulted  in  the  highest  difference  between 
MONO 7  and  NEC,  or  when  the  antenna  radius  was  0.02 a.  Table 
VII  shows  the  results  of  this  comparison  for  a  plate  spacing 
of  l . 4 A .  Once  again,  the  calculated  input  impedances  differ 
very  little  from  those  obtained  using  the  separate  gap 
expansion. 


Table  VII  COMPARISON  OF  M0N07,  M0N08  FOR  SELECTED  h 


Antenna 

M0N07 

M0N08 

%  Diff 

Length 

Rin 

Xin 

Rin 

Xin 

R 

X 

_ 

0 . 12  A 

13.35 

-92.19 

13.32 

-92.13 

0.22 

isa 

0 . 24  A 

46.08 

15.90 

46.08 

15.87 

0.00 

0 . 48  A 

220.12 

-120.91 

219.44 

-121.34 

0.31 

Ksa 
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As  noted  in  Chapter  4,  the  results  are  more  consistent 
with  the  solutions  obtained  from  NEC  when  Nx  equals  N2.  Figure 
25  illustrates  that  this  same  phenomenon  occurs  for  M0N08. 
Further  investigation  and  study  is  required  to  understand  the 
relationship  that  the  number  of  modes  in  Regions  I  and  II  have 
on  each  other.  Once  this  is  understood,  it  may  be  possible 


Comparison  of  the  Current  Distribution 


For  selected  N1«N2  respectively 


Figure  25  Current  Distribution  using  M0N08  as  Nx  and  Nz  Vary 


to  derive  an  expression  for  the  minimum  number  of  modes 
required  to  obtain  results  within  a  specified  accuracy  range. 
This  expression  could  then  be  used  to  allow  the  software  to 
dynamically  adjust  the  number  of  modes  required  to  obtain 
accurate  results  at  run  time. 


65 


In  all  the  above  cases,  the  results  were  compared  to  the 


solution  using  H0N07  where  Mx  was  chosen  as  21/ d.  However, 
Figures  21  and  22  (page  56)  Indicates  that  there  Is  an  optimum 
number  of  modes  for  the  gap  region  to  reduce  the  differences 
between  the  solutions  of  MONO  and  NEC.  Table  VIII  shows  that 
M0N08  returns  the  same  result  for  the  Input  Impedance  as  M0N07 
did  where  Mlf  Nx  and  Nz  were  equal.  While  this  result  may  not 
be  the  closest  to  NEC's  results,  the  differences  are  consis- 

Table  VIII  COMPARISON  OF  M0N08  AND  M0N07  AS  Mx  VARIES 


N1=N2 

M0N08 

R/X 

Ml=60 

M0N07  R/X 

Ml=50  Ml=40 

l 

Ml=30 

60 

43.31 

43.31 

43.02 

43.39 

42.44 

13.87 

13.87 

13.85 

13.95 

14.14 

50 

42.59 

42.81 

42.59 

43.12 

42.13 

13.42 

13.44 

13.42 

13.53 

13.72 

30 

40.70 

41.61 

41.66 

41.72 

40.70 

12.64 

12.32 

12.32 

12.39 

12.64 

tently  less  than  ten  percent. 

In  conclusion,  consistent  results  are  obtained  without 
the  separate  expansion  for  the  gap  region.  While  these 
results  may  not  offer  the  closest  possible  solution  to  that 
available  by  NEC,  the  small  difference  in  the  solution 
obtained  without  the  extra  terms  is  more  than  offset  by  the 
significant  reduction  in  required  memory  and  calculation  time. 
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When  the  relationship  between  the  physical  structure  and  the 
minimum  number  of  modes  that  are  required  to  obtain  the  best 
possible  solution  is  understood,  then  the  software  can  be 
modified  so  as  to  select  the  number  of  terms  required  at  run 
time.  Until  then,  M0N08  should  be  the  method  of  choice  as 
its  results  are  more  predictable  and  is  capable  of  returning 
the  result  within  acceptable  accuracy  limits  without  "trial 
and  error". 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

We  have  shown  that  one  can  calculate  the  current  dis¬ 
tribution  along  a  monopole  using  multiregional  cylindrical 
harmonic  expansions  that  are  consistent  with  the  results 
obtained  from  the  Numerical  Electromagnetics  Code  (NEC) . 
While  this  software  package  is  not  intended  as  a  replacement 
for  NEC,  it  is  the  beginning  step  in  the  development  of 
software  directed  towards  solid  antenna  structures  that  are 
rotated  about  a  symmetrical  axis  with  a  simple  and  flexible 
input  structure.  The  particular  software  that  has  been 
developed  has  the  following  characteristics: 

•  Simple  input  data  structure.  All  inputs  entered  direct¬ 
ly,  via  unformatted  ASCII  file,  or  redirected  using 
standard  DOS  redirection  codes  allowing  batch  processing 
of  multiple  runs. 

•  Very  flexible  for  a  given  monopole  antenna  structure. 

•  Can  be  readily  modified  to  allow  for  a  homogeneous 
dielectric  coating  on  the  surface  of  the  antenna.  This 
coating  could  be  used  to  modify  the  surface  currents  and 
the  far  field  radiation  pattern. 

As  currently  written  for  the  simple  monopole,  the  user  is 
able  to  modify  all  the  physical  parameters  of  the  structure 
including  driving  source  frequency.  With  some  modifications, 
including  the  incorporation  of  a  finite  element  algorithm  in 
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one  sub-region,  the  user  may  also  include  the  effect  of  a 
nonhomogeneous  dielectric  coating  so  that  the  far  field 
radiation  pattern  may  be  tailored  as  desired. 

One  difference  that  should  be  noted  between  the  techni¬ 
ques  used  by  NEC  and  MONO  is  the  location  used  to  calculate 
the  input  impedance.  Due  to  the  NEC's  segmentation  of  the 
entire  structure,  the  input  impedance  is  actually  calculated 
in  the  center  of  the  gap  region.  However,  MONO  calculates  the 
input  impedance  based  on  the  value  of  the  current  at  the  base 
of  the  antenna.  This  difference  may  lead  to  large  discrepan¬ 
cies  in  the  calculated  input  impedances  if  the  value  of  the 
current  is  changing  rapidly  at  the  end  of  the  gap  region,  as 
was  observed  for  monopoles  whose  height  approaches  a  half 
wavelength.  Other  differences  include: 

•  MONO  requires  longer  computation  time  for  larger  struc¬ 
tures  than  does  NEC. 

•  MONO  Requires  larger  available  memory. 

•  MONO  currently  has  limited  error  detection  and  correction 
within  the  subroutines. 

•  The  current  inability  to  calculate  a  number  of  internal 
parameters  requires  a  "trial  and  error"  approach  to 
determining  the  minimum  number  of  modes  required  for 
proper  convergence. 

•  User  is  required  to  make  extensive  modifications  to 
MONO's  code  for  various  types  of  antenna  structures. 

•  MONO  is  a  highly  specific,  concept  validation  platform 
vice  a  generally  applicable  software  package. 
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B.  RECOMMENDATIONS 

Continued  research  is  required  in  order  to  understand  the 
effect  that  the  apparent  relationship  between  Nx  and  N2  has  on 
the  accuracy  of  the  result.  Additionally,  there  appears  to 
be  an  optimal  number  of  terms  required  for  the  gap  region  that 
results  in  an  apparent  improvement  of  the  result  when  compared 
to  those  obtained  using  NEC.  These  unanswered  questions  point 
out  the  need  for  additional  analysis  of  the  relationship 
between  the  number  of  terms  required  to  represent  the  electric 
field  in  a  specific  region  and  the  overall  accuracy  of  the 
final  result.  When  this  is  accomplished,  modifications  should 
be  made  to  allow  the  software  to  calculate  the  optimal  number 
of  terms  required  for  each  expansion  that  would  yield  results 
in  the  shortest  time  possible  for  a  specified  accuracy  range. 

Another  modification  could  be  made  to  the  code  that  would 
reduce  the  amount  of  memory  required  to  run  the  program.  This 
can  be  accomplished  by  temporarily  storing  only  one  row  of  the 
Pnk  matrix  at  a  time,  vice  storing  the  entire  matrix.  This 
approach  could  reduce  the  amount  of  required  memory  by  as  much 
as  50  percent. 

To  better  assess  the  benefits  of  this  method  over  the  one 
used  in  other  computational  schemes,  experimental  data  should 
be  compared  to  the  computed  results  for  a  variety  of  struc¬ 
tures.  In  this  way,  a  basic  understanding  of  when  this  method 
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best  approximates  the  actual  performance  characteristics  of 
an  operational  antenna  can  be  obtained.  Once  these  areas  are 
better  understood,  a  collection  of  routines  for  various 
physical  structures  should  be  developed  using  the  basic 
technique  developed  in  Chapter  2 .  This  collection  could 
include  top  loaded  monopoles,  and  conical  and  other  structures 
that  feature  symmetry  about  a  axis  of  rotation. 
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APPENDIX  A 


SOFTWARE  FLOW  CHART  FOR  MONO  WITH  GAP  FIELD  EXPANSION 
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APPENDIX  B 


SOURCE  CODE  FOR  MONO  WITH  GAP  FIELD  EXPANSION 


PROGRAM  MONO 
C  (Revision  7.3) 

C 

C  Computing  Currents  on  a  Monopole  Between  Parallel 
C  Plates  using  Cylindrical  Hanaonic  Expansions  in  2  Regions. 

C  Links  with  SUBS4.06J.  Created  05/87  LT.  R.C.  Hurley 


C 

C  . . -  Variables  - 

C  FMHZ  *  Frequency  in  MHz 

C  h  *  Monopole  height  in  Meters 

C  a  s  Monopole  radius  in  Meters 

C  d  s  Distance  between  lower  groind  plane  and  feed  point 

C  in  Meters 

C  l  3  Distance  between  ground  ptanes  in  Meters 

C  q  3  height  fri>?  top  of  monopole  to  upper  ground  plane 

C  in  Meters 

C  j  3  imaginary  One 

C  pi  3  value  of  pi  (3.1415927) 

C  ICO  3  2*pi*F/c 

C  where  F  *  frequency  of  operation  in  Hz 

C  c  3  Speed  of  Light 

C  K2  *  ko**2 

C  N1.N2  *  Number  of  Coefficients  in  each  of  the  two  regions 

C  NBIG  >  Size  of  the  arrays  used  in  FACTOR  and  SOLVE 

C 

C  Unit  6  3  Screen  Output 

C  Unit  7  3  Main  Output  File  to  Disk  (Name  selected  by  User.) 

C  *****  variable  Definitions  ***** 

C 


INTEGER  k,m,M1,Mmax,mpn,n,ns,N1,N2,Nbig,NI,Nmax,P(61) 

REAL  a,ai,Am,An,Bl,c,Cl,CM,CP,CR,d,di,OeltaZ,Denom,eps,FMKZ,h,hi 
REAL  1(0:500, 1 :3), ICO, K2,  l, li, Lambda, Pi, q,Rin,S 
REAL  T(0:500,0:60),w,Xin,y,yhat 

COMPLEX  Amk(61,61),B(61),C1,C2,IZ(61),j,JO,41,Ja(0:500,1:4) 
COMPLEX  Pmk(0:60, 0:500) ,Su*,Sgni1  ,Y0,T1  ,Zin 
CHARACTER*25  Fname 
CHARACTERS  Bell 
CHARACTER*)  HOW 
CHARACTER*)  INPUT 
CHARACTERS2  INname 
COMMON/Plkl/I, j,N1,Pi,d,Nbig 
COMMON/Blk2/K2, l,Mmax,N2,q,a,T 
C 

C  ********************************* 

C  *  Main  Program  Code  * 

C  ********************************* 

C 

c  *********************** 

C  *  Input  Section  * 

C  *********************** 

c 

Bel l*CHAR(7) 
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Nbfg*61 

c-3.06+08 

WRITE(6,*)  'Do  you  wish  to  use  an  input  data  file  tY/N]?' 

READ(0, 100)  HOW 

IF  (HOW.EQ.'y'.OR.HOU.EO.'Y')  THEN 

l*ITE{6,*)  'Enter  input  data  file  name  with  extension  :' 

READ(0,100)  IN name 

0PEN(Unit*4,File*INname) 

READ(4,100)  INPUT 

IF  ((INPUT  .£Q.  'w')  .OR.  (INPUT  .EQ.  'W')  .OR.  (INPUT  .EO. 

♦  '*')  .OR.  (INPUT  .EO.  'M')>  THEN 

READ(4,*)  FMHZ,hi ,ai ,di , l i ,N1 ,N2,M1,NI 

ELSE 

URITE(*,*)  Bell 

HRlTE(6,*)'Error  detected  in  first  line  of  data  file  ' 

+,INn ame 

WR1TE(6,*)  'First  line  must  be  either  H  or  W  for  meter 
+  values  or  wavelength  values.' 

STOP 

END  IF 

ELSE 

2  WRITE(6,*>  'Enter  data  in  either  (MJeters  or  [W]avelength  [M 

+/W]?' 

READ(0,100)  INPUT 

IF  (INPUT. EQ.'m'. OR. INPUT. EQ.'N')  THEN 

WRIT£(6,*>  'Enter  Frequency  in  MHz:  ' 

READ(0,*)  FHHZ 

URITE(6,*)  'Enter  Monopole  Height,  h,  in  meters:  ' 
READ(0,*)  hi 

URITE(6,*)  'Enter  Monopole  Radius,  in  meters:  ' 
REAO(0,»)  ai 

WRITE16,*)  'Enter  Feed  Gap,  d,  in  meters:  ' 

READ(0,*)  di 

URITE(6,*)  'Enter  Upper  Plane  Height,  l  (meters):  ' 
READ(0,*)  li 

ELSEIF  (INPUT. EO.'w'. OR. INPUT. EO.'U')  THEN 

WRITE(6,+)  'Enter  Frequency  in  MHz:  ' 

READ(0,+)  FMHZ 

WRITE(6,*)  'Enter  Monopole  Height,  h,  as  factor  of  wav 

♦elength  :  ' 

REAO(0,*)  hi 

WRITE(6,*)  'Enter  Monopole  Radius,  as  factor  of  wavele 

♦ngth  :  ' 

READ(0,*)  ai 

URITE(6,*)  'Enter  Feed  Gap,  d,  as  factor  of  wavelength 

+  :  ' 

READ(0,*)  di 

URITE(6,*)  'Enter  Upper  Plane  Height,  l,  as  factor  of 
♦wavelength  :  ' 

READ{0,*)  li 

ELSE 

VRITE(*,*)  Bell 

WRITE(6,*)  'You  must  enter  either  an  H  or  U. 

♦Reenter..  ' 

GOTO  2 

END  IF 

WRITE(6,*)  'Enter  #  Region  1  Coeffs.  (N1  .LE.  60):  ' 

READ(0,*)  N1 

WRITE(6,*)  'Enter  §  Region  2  Coeffs.  (N2  .LE.  60): 

READ(0,*)  N2 

WRITE(6,*)  'Enter  «  Gap  Region  Coeffs. (Ml  .LE.  500):  ' 

READ(0,*)  Ml 

WR1TE(6,*)  'Enter  No.  Points  for  !(z)  including  ends:  ' 
READ(0,*)  NI 
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o  o  o  o  o  r>  ooooo 


END  IF 

IF  (M1.GT.500)  Then 
UR1TE(*,*)  Bell 

URITE(6,*)  'Ratio  l/d  too  large.  Adjust  either:' 
WRITE(6,*)  '  a)  l  or  d  such  that  2l/d  <  500' 

WRITER,*)  '  b)  Array  sites  of  Im,T,Ja' 

STOP 

END  IF 


C  * 

C  *  Calculate  the  Physical  Parameters 

C  *  in  both  Wavelenght  and  Meters 

C  * 


Larbda*c/(FHHZ*10**(6)) 

IF  (INPUT. EO.'u'.OR. INPUT. EQ.'U')  THEN 
h«hi*Lanbda 
a-ai*Lanbda 
d=di*Larnbda 
l=li*Lambda 

ELSE 

h=hi 

hi=h/Lambda 

a-ai 

ai*a/Lambda 

d=di 

di=d/Lamb.. 

1=1  i 

li-l/Lambda 

END  IF 

C  *********************************************** 

•  * 

*  Is  Upper  Plate  Height  an  * 

*  Even  Multiple  of  Wavelength?  * 

*  • 
*************************  ********************* 

IF  (H00(( 1*2) .Lambda). IT.  IE-05)  THEN 
WRITE(*,*)  Bell 
URITE(6, *)  '  ' 

WRITEC6,*)  'ERROR!  Upper  Plate  Height  is  almost  an  even  mul 
♦tiple  of  half  the  wavelength.' 

URITE(6,*)  'Ensure  l  is  NOT  within  0.00001  of  an  integer  mul 
+tiple  of  half  the  wavelength.' 

STOP 

END  IF 

•A********************************************* 

•  • 

*  Is  Upper  Plate  Height  an  * 

*  Even  Multiple  of  Antenna  Height?  * 

*  * 

*********************************************** 
c  E*(10*h)-(ANINT((10*h)/l )*l ) 

IF  «M0D(l,h)).LT. IE-05)  THEN 
WRITEl*,*)  Bell 
WRITE<6,*)  '  ' 

WRITE(6,*)  ’ERROR!  Upper  Plate  Height  is  almost  an  even  mul 
♦tiple  of  the  antenna  height.' 

WRlT".t,*)  ’Ensure  1  is  NOT  within  0.00001  of  an  integer  mul 
♦tiple  of  h  (in  wavelength).' 
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STOP 


END  IF 

URITE(6,*)  'Enter  output  file  name  and  extension:  ' 

READ(0, 100)  Fname 

OPEN  (Unit*7fFile*Fname) 

UR1TE(7,210)  FMH 2 
WRITE<7,220)  h,hi 
WRITE(7,230)  a,ai 
WR1TE(7,240)  d,di 
WRITE(7,250)  l,ti 

WRIT£<7,*)  'Region  I  Coeffs.  (Ml):  ',N1 
URITE(7,*)  'Region  II  Coeffs.  (N2):  ',N2 
URITE(7,*>  'No.  Gap  Coeffs.  (Ml):  ',M1 
URITE<7,*)  'No.  Points  for  I(z):  ',NI 

j*(0.0,1.0) 

Pi=3. 1415927 
K0=Pt*FMHZ/150.0 
K2=K0*K0 
q*l-h 

Nmex=MAX0(N1 , N2) 

Mmax=MAX0(Nmax,M1 ) 

URITE(6,*)  'loading  the  Moment  Matrices  .  ' 

WRITE(6,*)  '  ' 


c 

***** 

c 

* 

* 

c 

* 

Evaluating  Moment  Integrations 

* 

c 

* 

Over  Region  h  <  z  <  l. 

* 

c 

* 

* 

c 

*********************************************** 

T(0,0)=q 

DO  11  k=1,N2 

T(0,k)=0.0 

11 

CONTINUE 

DO  22  m=1,M max 
Am=m*Pi/l 

DO  22  n=0,N2 

An=n*Pi/q 

IF((ABS(Am-An)).GT.1.0E-3)  THEN 

Oenom*Am*Am-An*An 

S=SIN(Aro*h) 

T(m,n)=-Am*S/Denom 

ELSE 

mpn=ro*n 

T(m,n)=0.5*q*((-1)**mpn) 

END  IF 

22 

CONTINUE 

WR1TE<6,«)  ' 
URITE(6/)  '  ' 

Moment  Mel  ix  T  loaded.' 

C 

***** 

C 

* 

c 

* 

Loading  the  lie  Matrix  for  the 

* 

c 

* 

Largest  Nuifcer  of  Coefficients 

* 

c 

* 

* 

c 

***** 

00  33  m=0,Mmax 

IF(m.EO.O)  Then 
I(m,1)=l 
I(m,2)*q 
!<m,3)«d 

ELSE 
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o  o  o  o  o  r> 


t(m, 1)*l/2.0 
I(«,2)-q/2.0 
Aws^Pi/l 

I(m,3)*SIM(*m*d)/fm 

END  IF 

33  CONTINUE 

URITE<6,*>  '  Moment  Matrix  I  Loaded.' 

WRITE<6,*)  '  ' 

URITE(6,*>  'Loading  System  Matrix  . ' 

CALL  ALQAD( Ja,Pmk,Amk) 

WRITE(6,*)  'System  Matrix  Amk  Loaded.' 

UR1TE(6,*)  '  ' 

URITE(6,*)  'Loading  Driving  Vector  . ' 

CALL  SLaAD(FNHZ,Ml,Pmk,4a,B) 

URITEC6,*)  'Driving  Vector  6  Loaded.' 

WtITE<6,*>  '  ' 

WRITE(6,*)  'Solving  the  Linear  System  ..  Please  Wait  ....  ' 
Ns«N1+1 

CALL  FACTOR(Amk,P,Ns,Nbig) 

CALL  SOLVE(Amk,B,P,Ns,Nbig) 

WRITE(7,#)  'Listing  the  Expansion  Coefficients' 

UR1TE<6,*)  '  Listing  the  Expansion  Coefficients  . ' 

URITE(7,*)  'N  and  B(N>!  ' 

DO  44  N=0,N1 

WRITE<7,«)  N,B(N+1 ) 

44  CONTINUE 


1  Computing  Current  from  the  * 

1  Expansion  Coefficients  * 

i********************************************** 


URITE<6,»)  '  ' 

URITE(6,*)  'Computing  current  on  the  monopole . ' 

URITE(7,*)  'Listing  the  Current  on  the  Monopole' 

Del taZ*(h'd)/(Nl - 1 ) 

URITE(7,*)  'Position  Distance  I (z)  (A)' 

UR1TE17,*)  '  Number  (Meters)  Mag  Phase' 

eps*8.854E-12 

u=2.0E6*Pi*FMHZ 

yhat«w*eps 

y-d 

k=1 

55  Sun=(0.,0.) 

DO  66  n=0,N1 

Sum=$um+(B(n+1 )*Ja(n,3)*cos(n*Pi*y/l )) 

66  CONTINUE 

Sum1=(0. 0,0.0) 

DO  77  mxO.Ml 

Sun1=sum1+((I(m,3)*Ja(m,3)*C0S(m*pi*y/l>>/(Js(m,4)*I(m,1))> 
77  CONTINUE 

I Z(k)*2 . 0*Pi*a*(Sum- j*yhat*Su»1/d) 

CM*CABS(IZ(k)) 

CR*REAL(IZ(k)) 

CI«AIMAG(IZ(k» 

CP*180.*ATAN2(CI ,CR)/pi 
URITE(7,200)  k,y,CM,CP 
y»y*OeltaZ 
k«k+1 

IF<y.LE.(M,t.E-6))  GO  TO  55 
WRITE(6,*)  'Complete.' 
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**************************************1 


c 
c 

C  *  Computing  Input  Impedance 

C  * 

C  ************************************1 

Zin=(1.0,0.0)/IZ<1) 

Rin*REAL(Zin) 

Xin*AIMAG(Zin) 

UR1TE(7,*)  'Input  Resistance  *  ',Rin 
WRITE(7,*)  'Input  Reactance  ■  ',xin 


c  ***********************************  ************ 

c  *  * 

C  *  I/O  Fonsat  Statements  * 
C  *  * 
C  ********************************* ************** 


100  FORMAT (A) 

102  FORMAT(I5) 

103  FORMAT(£12.3) 

200  FORHAT(2X,l3,5X,F8.3,5X,1PE12.3,5X,0PF8.2,1P£12.3,5X,0PF8.2, 

♦1PE12.3,5X,0PF8.2) 

210  FORMAT ( '  Driving  Frequency,  ',9X,F7.3,'  MHz') 

220  FORMAT ( '  Monopole  Height,  ',9X,F7.3,'  meters  (',F5.3,'  Lambda)') 

230  FORMAT C '  Monopole  Radius,  ',9X,F7.3,'  meters  (',F5.3,'  Lanbda)') 

240  FORMAT ( '  Feed  Gap  Distance,  '.9X.F7.3,'  meters  (',F5.3,'  Lanbda)') 

250  FORMAT ( '  Upper  Plane  Height, ',9X,F7. 3, '  meters  C',F5.3,'  Lwbda)') 

STOP 


r 
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Subroutine  Source  Code  for  MONO  with  Gap  Field  Expansion 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  ALQAD(Ja,Pmk,Amk) 
(Revision  7.3) 


*  Variable  Definitions  * 

•  • 


l  *  Distance  between  ground  planes 
Radius  *  Radius  of  monopole 

q  *  Distance  between  top  of  monopole  and  t^per  plane 
d  *  Antenna  feed  gap 

K2  ■  K0**2  (From  main  program) 

T(m,n)  »  Moment  integration  result  as  function  of  m  and  n 

I(m,x)  ■  Moment  integration  result  of  order  x  given  m 

j  >  Imaginary  One 

z  a  intermediate  value  for  Bessel  Subroutine 

un,vn  «  Eigenvalues 

JO  =  Bessel  Function  of  the  First  Kind  of  Order  0 

J1  a  Bessel  Function  of  the  First  Kind  of  Order  1 

HO  a  Hanfcel  Function  of  the  Second  Kind  of  Order  0 

Hi  a  Hanlcel  Function  of  the  Second  Kind  of  Order  1 

SUM  a  simulation  over  n  for  the  Amk  matrix 

SI  a  individual  term  used  in  SUM 

m  a  column  position  of  matrix  Amk 

k  a  row  position  of  matrix  Amk 

Hi  a  Member  of  sifedi visons  in  Region  1 

M2  a  Member  of  subdivisons  in  Region  2 

Amk  «  Complex  Matrix  value  that  describes  the  system 

This  SUBROUTINE  loads  the  complex  A  matrix  that  describes  the  system 
of  size  N1  by  Nl.  Maximus  size  of  matrix  is  set  at  50  by  50.  This 
SUBROUTINE  calles  to  the  following  subroutines  found  in  SUBS. FOR: 

CSR  -  Complex  square  root  that  preserves  negative  imaginary 
portion. 

BES1  •  Computes  the  Bessel  Function  values  for  a  given  z. 

HAN1  -  Computes  the  Hankel  Function  values  for  a  given  z. 

Written  5/88  by  R.C.  Hurley  using  Microsoft  Fortran  4.01. 

INTEGER  k.m.Mmax.n.Nbig^l.NZ 

REAL  l, a,q,d,K2.T(0:500, 0:60), 1(0:500. 1:3), Pi, R1,R2 
COMPLEX  j,zv,zu,vn,un,J0n,J1n,H0n,H1n,SUMfS1 
COMPLEX  CSR,Afflk(Nbig,Nbig),Ja(0:500,1:4) 

COMPLEX  Pa*<0:60, 0:500) 
character*!  er 

COMMON/Blk1/l,j,N1,Pi,d.Nbig 
COMMON/Blk2/K2, l,Mmax,N2,q,a,T 
cr»char(13) 

This  section  loads  a  matrix  consisting  of  Bessel  fuctions  and 
other  values  required  to  load  the  Amk  matrix.  The  colunns  are 
filled  as  follows: 

Ja(n,1)  «  un*Jo(un*a)  in  numerator  of  Sun 
Ja(n,2)  «  Jl(un*a)  in  denominator  of  Sun 
Ja(n,3)  *  vn*Hl(vn*a)  used  for  all  locations 
Ja(n,4)  *  vn*vn*Ho(vn*a)  used  in  diagonal  elements 

URITE(6, 100)  Mmex 


80 


DO  11  n*0,Mmx 

R1«K2-((n*Pi/q)*(n*Pi/q>) 

R2»K2-((n*Pi/l)*(n*Pi/l)> 

un*CSR(R1) 

vn>CSR(R2) 

ru«un*« 

zv*vn*a 

Call  8ES1(zu,J0n,J1n) 

Call  HAN1(zv,H0nfH1n) 

Ja(n,1)  *  Un*J0n 
Ja(n,2)  *  Jin 
Ja(n,3)  *  Vn*H1n 
Ja(n,4)  *  Vn*Vn*H0n 

c  urite(7,*)  'n*',n,' ja(n,4}=', ja(n,4) 

11  Continue 

MRITE(6,*)  '  Completed.' 

WRITE(6,*)  '  ' 

WRITE(6,200)  n1,Hmax 
C 

C  Loading  the  Pmk  Matrix  (size  N1  x  Mmex) 

C 

DO  22  w=0,N1 

WRITE(6,300)  m,cr 
DO  22  k*0,Mmax 

Pmk(m,k)«(0. 0,0.0) 

DO  33  n=0,N2 

S1*(Ja(n,1)*T(m,n)*T(k,n))/(I(n,2)*Ja(n,2)) 

Pmk(m,k)=Pmk(m,k)+S1 
33  CONTINUE 

22  CONTINUE 

URITEC6,*)  '  Completed.  ' 

DO  44  m*0,N1 

DO  44  k=0,N1 

If(m.EQ.k)  Then 

Amk(ro*1,k+1)=Ja(k,3)*Pmk(m,m)-)B(k,4)*l(m,1) 

ELSE 

Amk(mH,k+1)=Ja(k,3)*Pmk(m,k) 

END  IF 


44 

CONTINUE 

100 

FORMAT ( ' 

Loading  Ja  matrix  of  size', 13,'  x  4 . ') 

200 

FORMAT (' 

Loading  Pmk  Matrix  of  size', 13,'  x',!3,'.... 

300 

FORMAT ( ' 
RETURN 

END 

Calculating  row',i3,a,\) 

C 

c 

SUBROUTINE  BL0AD(FMHZ,M1 ,Pmk, Ja,B) 


C 

C  *’ 

C  * 


C 

C 

* 

* 

Variable  Definitions  * 

* 

c 

c 

B(m) 

* 

Driving  Matrix 

c 

N1 

* 

Number  of  Coefficient  Expansions  in  Region  1 

c 

I3(m) 

x 

Trig.  Moment  Integrations 

c 

FMHZ 

at 

Driving  voltage  frequency  in  megahertz 

c 

epa 

x 

8.854E-12 

c 

j 

X 

Imaginary  One 

C 


Real  d, 1(0:500,' ,3),Pi,eps,w 
Integer  k,M,M1,N1,Nbig 

Complex  B(Nbig), j,Ja(0:500,1 :4),Pmk(0:60, 0:500), Sum 
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U  U  U  U  (_>(_)  CJ  u  u  u  u  u  u  o  u  o  uuu  «-  u  o  o 


CCMM0N/Blk1/IJ,Nl,P1,d,Nbig 
epa«8.854E-12 
w«2.0E6*Pi*FMHZ 
00  22  m-O.NI 

SusKO. 0,0.0) 

DO  33  k«0,M1 

SwmSw*<<I(k,3>*Ja(k,3>*Pmk(m.k))/<l<k,1)*Ja(k,4>)) 
33  CONTINUE 

B(«r*1  )«j*w*eps*Sum/d 
22  CONTINUE 
RETURN 
END 


FUNCTION  CSR(R) 

REAL  R 
COMPLEX  CSR 

Evaluating  SORT(R)  with  special  consideration  to  neg.  reals 
IF(R.GE.O.O)  then 
CSR*SQRT(R) 

else 

CSR*<0.,-1.)*SQRT(-R) 

end  if 
RETURN 
END 


A******************************************************************* 


SUBROUTINE  BES1(Z, JO, J1 ) 

Computing  Bessel  Functions  for  n=0,1  with 
Complex  Argunent,  Z.  Direct  Power  Series  Method  for 
CABS(Z)  .LE.  6  and  Kankel's  Asymptotic  Formula  for 
CABS(Z)  .GT.  6.  Written  11/5/87  by  N.A.  Morgan 

INTEGER  M.M2 

REAL  C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX  Z,Z2,Z3,Z«,JO,J1,AM,CL,PO,P1,QO,01,CO,C1,SO,S1 

Pi*3. 1415927 

P2*2.0/PI 

1F(CABS(Z).LE.6.0)  THEN 

Utilizing  the  Direct  Power  Series  Method 

G0=  1.781072 

Z2«0.5*Z 

CL”CLOG(G0*Z2) 

Computing  F(m)  *  n  I  and  P(m)  *  1  ♦  1/2  ♦  1/3  *  1/m 

F(1)«1.0 
P(1)*1.0 
DO  11  M*2,34 

F(M)»M*F(M-1) 

P(M)*P(M-1)+1.0/M 
1  CONTINUE 

Computing  Power  Series  Coefficients 


22 

C 

C 

c 


33 


C 

C 

C 


DO  22  M«1,34 

C(M>aOM/<F(M)*F(M)) 

DM*-DM 

CONTINUE 

Confuting  JO  and  Jl 

J0«(1.,0.) 

J1«(0.,0.) 

M=0 

M=M+1 

M2*2*M 

AMaC(M)*<Z2**M2> 

J0*J0+AM 

J1«J1-M*AM 

IF((CABS(AM).GT.1.0E-10).AND.<M.LT.34»  GO  TO  33 

J1«J1/Z2 

return 

ELSE 

Hankel'  Asymptotic  Formula  (Abram,  t  Stegun  p.  364) 

Z2=Z*Z 

Z3=Z*Z2 

Z4=Z*Z3 

P0=1 .0- . 07031 25/Z2+.1 121521 /Z4 
00*- .125/Z+.0732422/Z3 
Pl=1.0*.1171875/Z2-.1441956/Z4 
Q1*.375/Z- .10253906/Z3 
CO=CCOS(Z-.25*PI) 

S0*CSIN(Z-.25*PI) 

C1-CCOS(Z-.75*Pl) 

S1»CSIN(Z-.75*PI) 

AM=CSORT (P2/Z ) 

J0*AM*(P0*C0-O0*S0> 

J1=AM*(P1*C1 -Q1*S1 ) 

ENDIF 

RETURN 

END 


C 

C  a********************************************************** ********* 
C 


C 

c 

c 

c 

c 

c 


c 

c 

c 


SUBROUTINE  HAH1(Z,H0,H1) 

Computing  Hankel  Functions  For  n*0,1  with 
Complex  Argunent,  Z.  Direct  Power  Series  Method  for 
CABS(Z)  .IE.  5  and  Hankel's  Asymptotic  Formula  for 
CABS(Z)  .GT.  5.  Written  11/6/87  by  M.A.  Morgan 

INTEGER  M,M2 

REAL  C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX  Z,Z2,Z3,Z4,j0,J1,Y0,Y1,AMfCL,P0,P1,Q0fO1 

COMPLEX  E0fE1,X0,X1,H0,H1,j 

Pl=3. 1415927 

P2»2.0/PI 

i*<0.,1.) 

1F(CABS(Z).LE.5.0)  THEN 

Direct  Power  Series  Method 


GO*  1.78072 
Z2*0.5*Z 
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non  £  onn  w  nnniy  nnn-»  n  n  n 


Ci.*CL0C(GO*Z2) 


Confuting  F(m)  *  ■  I  and  P(m)  *  1  +  1/2  ♦  1/3  + _ ♦  1/m 

F(1)*1.0 
P<1)*1.0 
DO  11  N«2,34 

F(H)»H*F<H-1> 

P(H)«P(M-1)*1.0/H 
1  CONTINUE 


Computing  Power  Series  Coefficients 


ON*- 1.0 
00  22  N*1,34 

C(M)*0N/{F(H)*F(N)) 

DN*-DM 

?  CONTINUE 

Computing  JO  and  J1 


J0=(1.,0.) 

J1=(0.,0.) 

H=0 

5  M=N+1 

N2*2*M 

AM=C(H)*(Z2**H2) 

J0=J0+AM 

Jl*4l-H*AM 

If((CA8S(AN).GT.1.OE-1O).AN0.<M.LT.34))  GO  TO  33 
J1*J1/Z2 

Computing  TO  and  Y1 
M=0 

Y0*CL*J0 

Y1*Z2*Cl*J1-0.5*J0 
*  H=H+1 

H2«2*M 

AM*C(H )*P(N )*  <  Z2**M2 ) 

Y0*Y0-AN 

Yl*Yl+N*AM 

1F<(CABS(AM).GT.1.0E-10).AND.<M.LT.34»  GO  TO  44 

Y0«P2*Y0 

Yl*P2*Y1/Z2 

HOajO- j*Y0 

Hl*Jl-j*Yl 

RETURN 

ELSE 

Hankel'  Asymptotic  Formula  (Abram.  &  Stegun  p.  364 

22*2*2 

23*2*22 

Z4«Z*Z3 

P0*1.0-.0703125/Z2*.1121521/Z4 
00*- . 125/2+ .0732422/Z3 
P1«1.0+.1171875/Z2-. 1441956/24 
01O75/Z-. 10253906/23 
X0»(2-.25*PI> 

X1*(2-.75*PI) 

E0*CEXP(- j*X0) 

E1=CEXP(- j*X1) 

AM*CSQRT(P2/Z) 


84 


HO«AM*(PO- j*O0)*E0 
Hl«AM*(P1- J*Q1)*E1 


END  IF 

RETURN 

END 


C 

C 

C 


SUBROUTINE  FACTOR  (A,P,N,NMX) 

C  PERFORMING  LU-DECOMPOSITION  WITH  PIVOTING  ON  THE  A-MATRIX 
C  WRITTEN  BY  M.A.  MORGAN 
C 

INTEGER  I, J, JP1,K,P(NMX),PJ/PR/R,RM1/RP1 

REAL  DMAX.ELMAG 

COMPLEX  A(NMX,NMX), 0(120) 

DO  60  R*1,N 

DO  10  K«1,N 

10  D(K)»A(K,R> 

RH1*R-1 

IF(RMI.LT.I)  GO  TO  31 
DO  30  <1-1,RM1 
PJ“P(J) 

A(J,R)=0(PJ) 

D(PJ)*D(J) 

JP1=J*1 
DO  20  I=JP1 ,N 

20  D(I)=D(I)-A(I,J)*A(J,R) 

30  CONTINUE 

31  CONTINUE 
DMAX«CABS(D(R)) 

P(R)=R 

RP1=R*1 

IF(RPI.GT.N)  GO  TO  41 
DO  40  IxRPl.N 

ELNAG*CABS(D(I» 

IF(ELMAG.LT.DMAX)  GO  TO  40 
DMAX=ELMAG 

40  P(R)*I 

41  CONTINUE 
PR»P(R) 

A(R,R)=0(PR) 

D(PR)*0(R) 

IF(RPI.GT.N)  GO  TO  51 
DO  50  I-RP1.N 

50  A(t,R)=0(I)/A(R,R) 

51  CONTINUE 

60  CONTINUE 

RETURN 

END 

C 

C 

SUBROUTINE  SOLVE  (A,B,P.N,NMX) 

C  BACK  SUBSTITUTION  TO  INVERT  THE  LINEAR  SYSTEM 

C  WRITTEN  BY  M.A.  MORGAN 

C 

COMPLEX  A(NMX,NMX),B(NMX), Y(120),SUM 
INTEGER  I,IP1,J,K,P(NMX>,PP1 
DO  20  I«1,N 
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PPI*P(1) 

Y(1)«S(PPI) 

B(PPI)*B<1) 

IP1«I*1 

IF(IPI.GT.N)  GO  TO  11 
DO  10  J»IP1,N 

B( J)*B( J)-A(J, I )*Y(I ) 

10  CONTINUE 

11  CONTINUE 
20  CONTINUE 

DO  40  K*1,N 

I*M-K*1 

SUM*(0.,0.) 

IPI-If 1 

IF(IPI.GT.N)  GO  TO  31 
DO  30  J«IP1,N 

SUN-SUM»A(1,J)*B(J) 

30  CONTINUE 

31  CONTINUE 
B(I)*(Y(I)-SUM)/A<I,l) 

40  CONTINUE 

RETURN 
END 

C 

c  jut***************************************************************** 

c 

PROGRAM  PLOT 

C  Program  to  Plot  a  Solid  Line  Plot  with  Option  to  Overlay  a 
C  ♦+++  Line  Plot  for  Comparison.  Uses  PLOT2  Data  File  Format. 

C 

C  Original  MS-FORTRAN  Version  of  PLOT2  4/24/87  by  M.A.  Morgan. 

C  Mods:  5/30/87  XMIM;  11/3/87  RM-FORTRAN  2nd  Plot  Option. 

C 

C  INPUT  DATA  FORMAT 

C 

C  Unit  #3:  Solid  Line  Plot 

C 

C  TITLE1-  64  Space  Header 

C  N  -  0  Data  Points  15  (Integer*2) 

C  XM1N  -  Real  Min  X  value  E12.3 

C  XMAX  -  Real  Max  X  value  E12.3 

C  F1(N)  -  Input  Data  Array  E12.3 

C 

C  Unit  0 4:  ♦  Symbol  Plot 

C 

C  TITLE2-  64  Space  Header 

C  N  -  0  Data  Points  15  <lnteger*2) 

C  XMIN  •  Real  Min  X  value  E12.3 

C  XMAX  -  Real  Max  X  value  E12.3 

C  F2(N)  -  Input  Oata  Array  E12.3 

C 

CHARACTER* 1  YN,YN1,DUM,YN2, SYMBOL 
CHARACTER*4  LINE 
CHARACTER*?  CROSS 
CHARACTER*^  LTIT.CTIT 
CHARACTER*32  SCALE 
CHARACTERS  TITLE, FNAME 
REAL  X( 1025), FI ( 1025), F2( 1025) 

INTEGER*2  N , JROW, JCOL , I SYMBL , I  TYPE 
SCALE*'ie  the  Vertical  Scale  Multiplier' 

LINE*'---  ' 

CROSS*'  ♦+♦  ' 
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CALL  KEEP 

URITEC*,*)  'Press  RET  to  Exit  from  Screen  Plot' 

URITEC*,*)  'Printer  Hardcopy  ?  (Y/N):  ' 

READ<*,100)  YH 

URITEC*,*)  'Enter  Data  File  Name  for  -  Plot:  ' 

READC*, 100)  FNAME 
C  Reading  froai  Unit  #  3 
0PEN<3,FILE*FNAME) 

READ <3, 100)  TITLE 
READ(3,120)  N 
RE AD (3, 130)  XMIN 
READ<3,130)  XNAX 

URITEC*,*)  ' -  Plot  TITLE,  N,  XMIN,  XMAX:  ' 

URITEC*,*)  TITLE 
URITEC*,*)  N, XMIN, XMAX 

URITEC*,*)  'Enter  -  Plot  Caption  <16  char  max):  ' 

READC*, 100)  LTIT 
DX=<XMAX-XMIN)/<N- 1.0) 

FMIN«0.0 
FMAX=0.0 
DO  22  K*1,N 
X<K)«XMIN+<K-1.0)*DX 
READ(3,130)  F1(X) 

IF(F1(K).LT.FMIN)  FMIN=F1<IC) 

IF(F1(K).GT.FMAX)  FMAX=F1<K) 

22  F2<K)=0.0 

URITEC*,*)  'Comparison  Plot  (+***)  on  Same  Graph  ???  (Y/N>:  ' 

READ(*, 100)  YN1 

IF<CYN1.EQ.'N').0R.<YN1.EQ.'n'))  GO  TO  11 
URITEC*,*)  'Enter  Data  File  Name  for  ♦+++  plot:  ' 

READC*,100)  FNAME 
C  Reading  from  Unit  #4 

0PEN<4,FILE*FNAME) 

READC4.100)  TITLE 
READC4.120)  N 
READ (4, 130)  XMIN 
READC4.130)  XMAX 

URITEC*,*)  ’****  Plot  TITLE,  H,  XMIN,  XMAX:' 

URITEC*,*)  TITLE 
URITEC*,*)  N, XMIN, XMAX 

URITEC*,*)  'Enter  ++++  Plot  Caption  <16  char  max):  ' 
READ<*,100)  CTIT 
DO  27  K*1 ,N 
READC4, 130)  F2<K) 

IF<F20C).LT.FMIN)  FMIN=F2<K) 

IFCF2CK).GT.FMAX)  FMAX=F2<K) 

27  CONTINUE 
CLOSE  <4) 

11  CONTINUE 

CALL  QBEEP 

URITEC*,*)  'Enter  TITLE  for  Plot  <64  Char.  Max):' 

URITEC*,*)  '  ' 

READ  <*,100)  TITLE 
CLOSE  <3) 

IFCFMIN.GT.O.O)  FMIN=0.0 
IFCFMAX.LT, 0.f  FMAX*0.0 
C  Computing  Sea'  Factors  for  Vertical  Axis 

ABSNIN-ABSCFMIN) 

ABSMAX-ABSCFMAX) 

YMAX-AMAX1 CABSMIN, ABSMAX) 

NSCL-INTCLOGIOCYMAX)) 

IF  CYMAX.LT. 1.0)  NSCL-NSCL-1 

YSCL«10.**NSCL 

FMIN*FM!N/YSCL 
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FNAX*FNAX/YSCL 
ABSMIN«AB$NIN/YSCL 
ABSHAX*ABS*AX/YSCL 
DO  33  K*1,N 
FI (K)»F1  OC)/YSCL 
33  F2<K)«F2(IC)/YSCL 
YM1N-0.0 

IF(FMIH.EQ.O.O)  GO  TO  37 
35  YNIN«YNIN*1.0 

IF(ABSMIN.GT.YMIN)  GO  TO  35 
YNIN»YMN*FNIN/ABS»HN 
37  CONTINUE 
VMAX«0.0 

IF(FMAX.EQ.O.O)  GO  TO  41 
39  YMAX»YMAX+1.0 

IF(ABSHAX.GT.YNAX)  GO  TO  39 
YNAX*YNAX*  FMAX/ ABSMAX 
41  CONTINUE 

C  Calling  GRAFMATIC  Routines  and  Plotting  FI  Solid  Line  Graph 
ITVPE-1 
ISYNBL«-2 
CALL  0SH00EC6) 

CALL  OCMOV<0,21) 

URITEf*, 150)  YSCL 

CALL  QPTXT(32,SCALE,3,25,20) 

CALL  OPLOT(100,530,12,147,XHIN,XMAX,YH1N,YMAX,0.,0. ,0,1 . ,1 .5) 
CALL  QSETUP(0,3, ISYHBL.3) 

XMA JOR*XMAX/5 . 0 

CALL  QXAXIS(XMN,XMAX,XNAJOR,1,1,2) 

YMAJOR-I.O 

CALL  OYAXIS(YNIN,YHAX,YHAJOR,1 ,1 ,2) 

CALL  0PTXT<64,TITLE,3, 14,24) 

CALL  QPTXT(4,LINE,3,14,22) 

CALL  QPTXT(16,LTIT,3,18,22) 

JRCU-<ABS<YMIN)/<ABS<Y«AX)+ABSCYNIN)))*135 

JCOL*70-430*XHIN/(XHAX-XNIN) 

CALL  QGTXT(3, '0.0' ,3, JCOL, JROW.O) 

CALL  QTABL(1TYPE,N,X,F1) 

IF(( YN1. EO.'N'). OR. (YN1 .EO.'n' ))  GO  TO  43 
C  Plotting  F2  Grapn  u?ing  +  Symbol 

CALL  OPTXT(7,CROSS,3,34,22) 

CALL  0PTXT(16,CTIT,3,41 ,22) 

ISYNBL«43 
I TYPE *0 

CALL  QSETUP(0,3,ISvnBL,3) 

CALL  OTABL(ITYPE,N,X,F2) 

43  CONTINUE 
READ(*,100)  OUM 

lF«YN.EO.'N').OR.(YN. EO.'n'))  GO  TO  44 
CALL  QPSCRN 

44  CONTINUE 

CALL  QSM0DE(2) 

URITEt*,*)  'End  of  PLOT  Progr*«' 

100  FORHAT(A) 

120  FORMAT(I5) 

130  FORMATCE12.3) 

150  FORMAT(13X,1PE11 .3) 

STOP 

END 
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SOFTWARE  FLOW  CHART  FOR  MONO  WITHOUT  GAP  FIELD  EXPANSION 
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APPENDIX  D 


SOURCE  CODE  FOR  MONO  WITHOUT  GAP  FIELD  EXPANSION 

PROGRAM  MONO 
C  (Revision  8.1) 

C 

C  Mod  of  Mono  7.2,  not  using  the  expansion  for  the  Gap  Region 
C  Computing  Currents  on  a  Monopole  Between  Parallel 
C  Plates  using  Cylindrical  Hanaonic  Expansions  in  2  Regions. 

C  Links  with  SU8S4.0BJ.  Created  05/88  LT.  R.C.  Hurley 


C  Modified  10/88  N.  A.  Morgan 

C 

C  .  Variables  - 

C  FMHZ  *  Frequency  in  MHz 

C  h  *  Monopole  height  in  Meters 

C  a  *  Monopole  radius  in  Meters 

C  d  *  Distance  between  lower  ground  plane  and  feed  point 

C  in  Meters 

Cl  *  Distance  between  ground  planes  in  Meters 

C  q  s  height  from  top  of  monopole  to  upper  gromd  plane 

C  in  Meters 

C  j  *  Imaginary  One 

C  Pi  *  Value  of  pi  (3. 1415927) 

C  KO  =  2*pi*F/c 

C  where  F  *  frequency  of  operation  in  Hz 

C  c  =  Speed  of  Light 

C  K2  «  K0**2 

C  HI, M2  ■  Hinder  of  Coefficients  in  each  of  the  two  regions 

C  MBIG  =  Size  of  the  arrays  used  in  FACTOR  and  SOLVE 

C  Unit  0  *  Standard  Input  Device 

C  Unit  6  *  Screen  Output 

C  Unit  7  *  Main  Output  File  to  Disk  (Name  selected  by  User.) 

C  *****  variable  Definitions  ***** 

C 

$Large:Amk 


IHTEGER  k,mtfflpn,nlN1,N2,NbiglNI,Nmax,Ns,N2,P(61) 

REAL  a,ai,Am,An,c,CI,CM,CP,CR,d,di,DeltaZ,Denom,Dz,EzM,eps,FMHZ 

REAL  h, hi, 1(0:60,1:3), K0,K2, l, l i, Lambda, Pi, q,Rin,S 

REAL  T(0:60,0:60),w,Xin,y,yhat,Z,Zmin,Zmax 

COMPLEX  Amk(61,61),B(61),C1,C2,IZ(61), j, J0,J1, Ja(0:60,1:4) 

COMPLEX  Pfflk(0:60, 0:60), Sun, Y0,Y1,Zin 

CHARACTER*25  Frame 

CHARACTER*1  Bell 

CHARACTERS  HOW 

CHARACTERS  INPUT 

CHARACTER*)  Eok 

CHARACTERS2  INnane 

COMMON/Blkl/I , j ,N1 ,Pi ,d,Mbi  g 

C0MM0N/Blk2/K2,l,Nmax,N2,q,a,T 


Main  Program  Code 
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Input  Section 


8elt«CHAR<7> 

Nbig-61 

c*3.0E+08 

WRITEC6,*)  'Do  you  wish  to  use  an  input  data  file  (Y/N)?' 
REMX0,100)  HOU 

IF  (HGW.EQ.'y'.OR.NOU.EQ.'Y')  THEN 

URITE(6,*)  'Enter  input  data  file  name  with  extension 
REMHO,  100)  INnaiae 
0PEN(Unit*4,Fi  le-IMnnw) 

REMX4,100)  INPUT 

IF  ((INPUT  .EO.  'M')  .OR.  (INPUT  .EO.  'M')  .OR.  (INPUT  .EO. 

♦  '•')  .OR.  (INPUT  .EO.  'M')>  THEN 

REMX4,*)  FMHZ,hi ,ai ,di , li ,N1 ,H2,NI 

ELSE 

WRITE(*,*)  Bell 

WRITE(6,*)'Error  detected  in  first  line  of  date  file  ' 

+,INname 

URITE(6,*)  'First  line  must  be  either  H  or  W  for  meter 

♦  values  or  wavelength  values.' 

STOP 

ENDIF 

ELSE 

2  WRITE(6,*)  'Enter  data  in  either  Meters  or  Wavelength  [M 

*/\nv 

REMX0.100)  INPUT 

IF  ( INPUT. EQ. 'at' .OR. INPUT. EO. 'M' )  THEN 

URITE(6,*)  'Enter  Frequency  in  MHz:  ' 

REMHO,*)  FMHZ 

WRITE(6,*)  'Enter  Honopole  Height,  h,  in  aieters:  ' 
REMHO,*)  hi 

URITE(6,*)  'Enter  Monopole  Radius,  in  aieters:  ' 
REMHO,*)  ai 

URITE(6,*)  'Enter  Feed  Gap,  d,  in  aieters:  ' 

REM)(0,*)  di 

URITE(6,*)  'Enter  Upper  Plane  Height,  l  (meters):  ' 
REMHO,*)  li 

ELSE IF  (INPUT. EO.'w'. OR. INPUT. EO.'U')  THEN 

WRITE(6,*)  'Enter  Frequency  in  MHz:  ' 

REMHO,*)  FMHZ 

MRITE(6,*)  'Enter  Monopole  Height,  h,  as  factor  of  wav 

♦elength  :  ' 

REMHO,*)  hi 

UR1TE(6,*)  'Enter  Monopole  Radius,  as  factor  of  wavele 

+ngth  :  ' 

REMHO,*)  ai 

URITE(6,*)  'Enter  Feed  Gap,  d,  as  factor  of  wavelength 

+  :  ' 

REMHO,*)  di 

UR!TE(6,*)  'Enter  Upper  Plane  Height,  l,  as  factor  of 
♦wavelength  :  ' 

REMHO,*)  li 

ELSE 

URITE(*,*)  Bell 

URITE(6,*)  'You  must  enter  either  an  M  or  U. 

♦Reenter..  ' 

GOTO  2 

ENDIF 

URITE{6,*)  'Enter  #  Region  1  Coeffs.  (Ml  .LE.  60):  ' 

REMHO,*)  N1 
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WRITE(6,*)  'Enter  §  Region  2  Coeffe.  (N2  .IE.  60):  ' 

READ(0,*>  N2 

WRITE(6f»)  'Enter  No.  Point*  for  I(i)  including  ends:  ' 
REAO(0,*)  MI 

ENDIF 


Calculate  the  Physical  Paraawters 
in  both  Uavelenght  and  Heters 


Lambda«c/(FHHZ*10**(6)) 

IF  (INPUT. EO.'w' .OR. INPUT.EO.'W')  THEN 
h*hi*Laafcda 
a«ai ‘Lambda 
d*di*Laafcda 
l=li*Lanbda 

ELSE 

h'hi 

hi^h/Lambda 

a=ei 

a i =8/ Lambda 
d*di 

di*d/Lambda 

1*11 

IM/Lambda 

ENDIF 


i*********************************************- 


Is  Upper  Plate  Height  an  < 

Even  Multiple  of  Wavelength?  1 

>  1 

r*********************************************i 


IF  (M00((l*2). Lambda). LT.  IE-05)  THEN 
URITE(*,*)  Bell 
WRITE(6,*)  '  ' 

WRITEC6,*)  'ERROR!  Upper  Plate  Height  fs  almost  an  even  mul 
♦tiple  of  half  the  wavelength.' 

URITE(6,*)  'Ensure  l  is  NOT  within  0.00001  of  an  integer  mul 
♦tiple  of  half  the  wavelength.' 

STOP 

ENDIF 


C 

C 

c 

c 

c 

c 

c 


Is  Upper  Plate  Height  an 
Even  Multiple  of  Antenna  Height? 


IF  <(M0D(l,h)).LT. IE-05)  THEN 

WRITE(V)  Bell 

URITE<6,*)  '  ' 

UR1TE(6,*)  'ERROR!  Upper  Plate  Height  is  almost  an  even  mul 
♦tiple  of  the  antenna  height.' 

UR1TE(6,*)  'Ensure  l  is  NOT  within  0.00001  of  an  integer  mul 
♦tiple  of  h  (in  wavelength).' 

STOP 

Em  i  f 

WRITE16,*)  'Enter  output  file  name  and  extension:  ' 
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READ(O.IOO)  Fname 
OKU  (Uni  t»7,  F i  l e« Frame) 

UR1TE(7,210)  FHHZ 
URITE(7,220>  h,hi 
WRITE(7,230)  a.ai 
WRJTE(7,240)  d,di 
WRITE (7, 250)  1,11 

WRITE{7,*)  'Region  l  Coaffa.  (Ml):  '.Ml 
WRITE(7,*)  'Region  II  Coaffa.  (M2):  ',N2 
WR1TE(7,*)  'Mo.  Points  for  I(i):  '.MI 

j»(0.0,1.0) 

Pi«3. 1415927 

K0*Pi*FMHZ/150.0 

K2*K0*K0 

epe-6.854E-12 

w*2.066*Pi*FMHZ 

yhat*w*epe 

<ri-h 

Mmex*MAX0(N1,M2) 

WRITE(6,*)  'Loading  the  Moment  Matrices  , 
WRITE(6,*)  '  ' 


C 

C 

C 

C 

C 

C 


****************** 


Evaluating  Moo.tnt  Integrations 
Over  Region  h  <  z  <  l. 


T(0,0)«q 
00  11  k*1.N2 

T(0,k)«0.0 
11  CONTINUE 

00  22  m*1,Nmax 
RwaPPi/I 

00  22  n*0,N2 

An*n*P1/q 

I F ( ( A8S( Am- An) ) . GT . 1 . OE -3)  THEN 
Denc«mAm*Am-An*An 
S=slM(Am*h) 

T(m, n)=-Am*S/Denom 

ELSE 

«pn*n*n 

T(a,n)»0.5*q*((-1)**irpn) 

END  IF 

22  CONTINUE 

WIITE(6,*)  '  Moment  Matrix  T  Loaded.' 

URITE(6,*)  '  ' 

************************** ********** 


Loading  the  Im  Matrix  for  the 
Largest  Member  of  Coefficients 


DO  33  m*0,Mmax 

IF(a.EQ.O)  Then 
I(M,1)*l 
I(a,2)*q 
I(m,3)*d 

ELSE 

l(m,1)*l/2.0 

l(m,2)*q/2.0 
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Amn*Pi/l 

1  (•,3)»SIN(Ajn*d)/AJH 

END  IF 

33  CONTINUE 

URITE(6,*)  '  Moment  Matrix  1  Loaded.' 

WRITE<6,*)  '  ' 

WRITE(6,*>  'Loading  System  Matrix  . ' 

CALL  ALOAD(Ja,Pnk,Aak) 

WRITE(6,*>  'System  Matrix  Amk  Loaded.' 

WRITE<6,*)  '  ' 

l*ITE<6,*)  'Loading  Driving  Vector  . ' 

DO  40  m>0,N1 

B(»*1)ej*yhat*I(a,3)/d 
40  CONTINUE 

URITE(6,*)  'Driving  Vector  B  Loaded.' 

URITE(6,*)  '  ' 

WRITE(6,*)  'Solving  the  Linear  System  ..  Please  Wait  ....  ' 
Ns«N1+1 

CALL  FACT0R(Amk,P,N8,Nbig) 

CALL  SOLVE(Amk,B,P,Ns.Nbig> 

URITE(7,*)  'Listing  the  Expansion  Coefficients' 

URITE(6,*)  '  Listing  the  Expansion  Coefficients  . ' 

URITE(7,*>  'N  and  B(N):  ' 

DO  44  N=0,N1 

URITE(7,*)  N,B(N+1) 

44  CONTINUE 


*********************************************** 

•  , 

*  Computing  Current  from  the  * 

*  Expansion  Coefficients  * 

*  * 
*********************************************** 

URITE16,*)  '  ' 

WRITE(6,*)  'Computing  current  on  the  monopole . ' 

URITE(7,*)  'Listing  the  Current  on  the  Monopole' 
DeltaZ*(h*d)/(NI-1) 

URITE(7,*)  'Position  Distance  I(z>  (A)' 

UR1TE(7,*)  '  Number  (Meters)  Mag  Phase' 

V*d 

k*1 

j  Sua=(0.,0.) 

DO  66  n*0,N1 

Sua=Sum+(B(n*1 )*Ja(n,3)*cos(n*Pi*y/l ) > 

S  CONTINUE 

IZ(k)*2.0*Pi*a*Sun 

CM-CABS(IZ(k» 

CR«REAL(IZ(k)> 

Cl*AIMAG(IZ(k)) 

CP*180.*ATAN2(CI , CR )/pi 
URITE(7,200)  k.y.CM.CP 
y»yK)eltaZ 
k«k+1 

IF(y.LE.<h*1.E-6)>  GO  TO  55 
WRITE(6,*)  'Complete.' 

******* **************************************** 

*  * 

*  Computing  Input  Impedance  * 

*********************************************** 

Zin*(1.0,0.0)/lZ(1) 
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Rln-REAL(Zln) 

Xin*AlNAG(Zin) 

URITE(7,*)  'Input  Resistance  *  ',Rin 
WR1TE(7,*>  'Input  Reactance  «  ',Xin 


Checking  Ez(a,z)  of  solution  at  200  points  from  0  to  L 
UR1TE(6,*)  'Do  you  want  to  compute  the  Electric  Field  (Y/N)?  ' 
READ(O.IOO)  Eok 

IF  ((Eok  .EG.  'Y')  .OR.  (Eok  .EG.  'y'))  THEM 
NZ-200 
ZMIN-0.0 
ZMAX-L 

DZ«ZMAX/(NZ-1.0) 

OPEN(3,FlLE*'EZ.DAT') 

URITE(3,100)  'EZ  vs.  Z' 

WRITE(3,102)  NZ 
URITE(3,103)  ZNIN 
WRITE(3, 103)  ZMAX 
DO  99  K*1,NZ 

Z*(K-1 . )*DZ 
C1*(0.,0.) 

00  88  N=0,N1 

C1=C1+B(N*1 )*Ja(N,4>*COS(N*Pi*z/l ) 

68  CONTINUE 

C1«C1/(]*yhat) 

EzM*CABS(C1) 

WRITE(3, 103)  EzM 

99  CONTINUE 

URITE(6,*>  '  Electric  Field  data  stored  in  file  EZ.DAT.' 

ELSE 

URITE(6,*>  '  Completed.' 

END  IF 
C 

^  *********************************************** 

C  *  * 

C  *  I/O  Format  Statements  * 

C  *  * 

C  *********************************************** 

100  FORMAT(A) 

102  FORNAT(IS) 

103  FORHAT(F12.3) 

200  FORMAT(2X,l3,5X,F8.3,5X,1PE12.3,5X,0PF8.2,1PE12.3,5X,0PF8.2, 

♦1PE12.3,5X,0PF8.2) 

210  FORMAT ( '  Driving  Frequency,  ',9X,F7.3,'  MHz') 

220  FORMATS  Monopole  Height,  '.9X.F7.3,'  meters  (',F5.3,»  Lambda)') 

230  FORHAT('  Nonopole  Radius,  '.9X.F7.3,'  meters  (',F5.3,'  Lambda)') 

240  F0RMAT('  Feed  Gap  Distance,  '.9X.F7.3,'  meters  (',F5.3,'  Lmibda)') 

250  FORMATS  Upper  Plane  Height, '.9X.F7.3,'  meters  (',F5.3,'  Lambda)') 

STOP 
END 
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Subroutine  Sorce  Code  for  MONO  without  Gap  Field  Expansion 

SUBROUTINE  ALQAD< Ja.Pmk.Ank) 

(Revision  B.1) 


Variable  Definitions 


C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

(LARGE: Amk 

INTEGER  k.m,n,Nbig.Nmax,Nl,N2 
REAL  l,a,q,d,K2,T(0:60,0:60),I<0:60,1:3>,Pi,R1,R2 
COMPLEX  j.zv.zu.vn.un.JOn.JIn.HOn.HIn.SUM.Sl 
COMPLEX  CSR.AmkfNbig.NbigMefOidO.I.'A) 

COMPLEX  Pmk(0:60,0:60) 

CHARACTERS  cr 
CCMN0N/Blk1/l,j,N1,Pifd,Nbig 
COMMON/ B l k2/K2, I ,Nmax,N2,q,a,T 
cr«CHAR<13) 


l  *  Distance  between  ground  planes 
Radius  ■  Radius  of  monopole 

q  «  Distance  between  top  of  monopole  and  upper  plane 
d  b  Antenna  faed  gap 

K2  a  xo**2  (Froai  Min  program) 

T(m,n)  »  Moment  integration  result  as  function  of  m  and  n 

I(m,x)  b  Moment  integration  result  of  order  x  given  m 

Imaginary  One 

Intermediate  value  for  Bessel  Subroutine 
Eigenvalues 

Bessel  Function  of  the  First  Kind  of  Order  0 
Bessel  Function  of  the  First  Kind  of  Order  1 
Hankel  Function  of  the  Second  Kind  of  Order  0 
Hankel  Function  of  the  Second  Kind  of  Order  1 
Simulation  over  n  for  the  Amfc  matrix 
Individual  term  used  in  SUM 
Column  position  of  matrix  Amk 
Row  position  of  matrix  Amk 
Number  of  subdivisons  in  Region  1 
Number  of  subdivisons  in  Region  2 
Complex  Matrix  value  that  describes  the  system 


J 

z 

un,vn 

JO 

J1 

HO 

HI 

SUM 

SI 

m 

k 

HI 

N2 

Amk 


This  SUBROUTINE  loads  the  conplex  A  matrix  that  describes  the  system 
of  size  Nl  by  N1.  Maximum  size  of  matrix  is  set  at  50  by  50.  This 
SUBROUTINE  calles  to  the  following  subroutines  found  in  SUBS. FOR: 

-  Complex  square  root  that  preserves  negative  iMginary 
portion. 

•  Computes  the  Bessel  Function  values  for  a  given  z. 

*  Computes  the  value  of  the  moment  integratation  T 

-  Computes  the  value  of  the  moment  integratation  I 


CSR 

BES1 

Tmint 

Imint 


Written  5/88  by  R.C.  Hurley  using  Microsoft  Fortran  4.01. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


This  section  loads  a  matrix  consisting  of  Bessel  fuctions  and 
other  values  required  to  load  the  Amk  matrix.  The  colums  are 
filled  as  follows: 

Ja(n,1)  ■  un*Jo(un*a)  in  numerator  of  Sua 
Ja(n,2)  >  Jl(un*a)  in  denominator  of  Sim 
Ja(n,3)  «  vn*H1(vn*a)  used  for  all  locations 
Ja(n,4)  «  vn*vn*Ho(vn*a)  used  in  diagonal  elements 
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UR1TE(6,100)  Mkx 
DO  11  n>0,NMX 

R1"K2*(<n*Pi/q)*(n*Pi/q)) 

R2»K2"«n*Pi/l)*<n*Pi/l)) 

un-C$R(R1> 

vn-CSR(R2>  , 

zu«un*a  & 

zv*vn*a 

Call  BES1(zu,J0n,J1n) 

Call  HAN1(zv,H0n,H1n) 

Ja(n,1)  ■  Un*J0n 
Ja(n,2)  *  Jin 
Ja(n,3)  ■  Vn*H1n 
Ja(n,4)  ■  Vn*Vn*H0n 
11  Continue 

URITE(6,*)  *  Completed.' 

URITE16,*)  '  ' 

URITE(6,200)  NI.Naau 
C 

C  Loading  the  Pa*  Matrix  (size  N1  x  Maax) 


DO  22  w*0,N1 

URITE(6,300)  *,cr 
DO  22  k*0,Na»x 

Pisk(ai,k)*<0. 0,0.0} 

DO  33  n*0,N2 

SWJa(n,1)*T(ai,n>*T(k,n»/(I(n,2>*Ja(n,2» 
Pa*(a,k)»Piak(»,k)+S1 
33  CONTINUE 

22  CONTINUE 

URITE16,*)  '  Coapleted.  ' 

DO  44  a»0,N1 

DO  44  k«0,N1 

If(M.EO.k)  Then 

Ank(aiM,k*1>*Ja(k,3)*PiBk<ai,ai)-Ja(k,4)*I(B,1) 

ELSE 


Aaklanl ,  k+1  )=Ja(k,3)*Piak<Bi,k) 

END  IF 

44  CONTINUE 

100  FORMAT ( '  Loading  Ja  aiatrix  of  size', 13,'  x  4 . ') 

200  FORHAT('  Loading  Pa*  Matrix  of  size', 13,'  x',I3,' . ') 

300  F0RMAT('  Calculating  row',i3,a,\) 

RETURN 

END 


C 

c  ************************************************************* ******* 

c 


FUNCTION  CSR(R) 

REAL  R 
COMPLEX  CSR 

C  Evaluating  SQRT(R)  with  special  consideration  to  neg.  reals 
IF(R.GE.O.O)  THEN 
CSR*SQRT(R) 

ELSE 

CSR«(0.,-1.)*SORT(-R) 

END  IF 
RETURN 
END 
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c 

c 

SUBROUTINE  BESKZ,J0,4l> 

C 

C  Computing  Bessel  Functions  for  n*0,1  with 
C  Complex  Arguaent,  Z.  Direct  Power  Series  Method  for 
C  CABS(Z)  .IE.  6  and  Henkel's  Asymptotic  Fonaula  for 
C  CABS(Z)  .CT.  6.  Written  11/5/87  by  N.A.  Morgen 
C 

INTEGER  M,M2 

REAL  C(34>,DMfF(34>,G0.P<3*),PifP2 

COMPLEX  Z,Z2,Z3,Z4,J0,  Jl,AM,CL,P0,P1,O0,Ql,C0,C1 ,S0,S1 

Pi «3. 1415927 

P2*2.0/PI 

1F(CABS(Z) .LE.6.0)  THEN 
C 

C  Utilizing  the  Direct  Power  Series  Method 
C 

GO*  1.781072 

Z2«0.5*Z 

CL*CLOG(GO*Z2) 

C 

C  Computing  F(m)  *  m  I  and  P(m)  «  1  ♦  1/2  +  1/3  ♦  ....♦  1/m 
C 

F(1)«1.0 
P(1)*1.0 
DO  11  M»2,34 

F<M)*N*F(M-1> 

P(M)»P(M- 1 )+1 . 0/M 
11  CONTINUE 

C 

C  Computing  Power  Series  Coefficients 

C 

DM«-1.0 
DO  22  M-1,34 

C(M)«OM/{F(M)*F(M)) 

DM*-DM 

22  CONTINUE 

C 

C  Computing  JO  and  J1 

C 

J0*(1.,0.) 

J1«(0.,0.) 

M-0 

33  M«M+1 

M2*2*M 

AMIC<M)*(Z2**M2) 

JO-JO+AM 

J1»J1-M*AM 

IF((CABS(AM).GT. 1 .0E-10).AM).<N.LT .34))  GO  TO  33 

J1-J1/Z2 

return 

ELSE 

C 

C  Henkel'  Asymptotic  Formula  (Abram.  <  Stegun  p.  364) 

C 

Z2«Z*Z 

23«Z*Z2 

Z4«Z*Z3 

P0»1.0-.0703125/Z2+.1121521/Z4 
Q0»- . 125/Z+.0732422/Z3 
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U  (J  (J  U  U  U  U  u  u  OUO  U  U  U  r-  U  U  U  OIU 


P1«1 .0O171875/Z2-  .1441956/Z4 
01* ,375/Z- . 10253906/Z3 
CO»CCOS(Z- .25*PI > 
S0*CSIN(Z-.25*Pl) 
C1«CC0S<Z*.75*PI) 

S1*C$IM(7  .75*Pl) 

AN-CS0RT(P2/Z) 

JO*AM*(PO*CO-GO*SO> 

r  am»<pi*ci-oi*si) 

END  IF 

RETURN 

END 


SUBROUTINE  HAN1(Z,H0,H1) 

Computing  Hank* l  Function*  for  r*«0,1  with 
Complex  Argument,  Z.  Direct  Power  Series  Method  for 
CABS(Z)  .IE.  5  and  Nankel's  Asymptotic  Formula  for 
CABS(Z)  .GT.  5.  Written  11/6/87  by  M.A.  Morgan 

INTEGER  M,H2 

REAL  C<34),DMfF(34),G0,P(34),Pi,P2 

COMPLEX  Z,Z2.Z3,Z4,JO,J1,YO.Y1,AM,CL,PO,P1,QO,01 

COMPLEX  E0,E1,X0,X1,N0,H1,j 

Pl*3. 1415927 

P2*2.0/PI 

j-<0.,1.) 

IF(CABS(Z).LE.5.0)  THEN 

Direct  Power  Series  Method 

GO*  1.78072 
Z2*0.5*Z 
CL*CLOG( G0*Z2 ) 

Computing  F(m)  *  m  I  and  P(m)  *  1  ♦  1/2  ♦  1/3  ♦  _ ♦  1/m 

F(1)*1.0 
P<1)*1.0 
DO  11  M-2,34 

F(M)*M*F(M-1) 

P(M)*P(M-1)*1.0/M 
1  CONTINUE 

Computing  Power  Series  Coefficients 


DM— 1.0 
DO  22  M*1,34 

C(M)«0M/(F(M)*F(M)) 
DM— DM 

2  CONTINUE 

C  Computing  JO  end  J1 
C 

J0*(1.,0.) 

J1*(0.,0.) 

M*0 

33  M-M+1 

M2«2*M 

AM*C(M)*<Z2**M2) 


100 


o  r>  non  oor>  t  o  o  o 


JO-JO+AM 

J1«J1-M*AM 

IF((CABS(AM).GT.1.0€-10).AND.(M.LT.34))  GO  TO  33 
41«J1/Z2 

Computing  YO  and  Y1 


M=0 

YO«CL*JO 

Y1*Z2*CL*J1*0.5*J0 

IWH1 

M2«2*N 

AM*C(M)*P(M)*(Z2**M2) 

Y0*Y0-AM 

Y1*Y1+M*AM 

IF((CABS(AM).GT.1 .0E-10).AND.(M.LT.34))  GO  TO  44 

Y0«P2*Y0 

Y1«P2*Y1/Z2 

HO*JO- j*Y0 

H1*J1-j*Y1 

RETURN 

ELSE 

Hankel'  Asymptotic  Formula  (Abram.  &  Stegun  p.  364 

Z2=Z*Z 

Z3=Z*Z2 

Z4=Z*Z3 

P0=1 .0- . 07031 25/Z2+.1 121521/24 
00= - . 1 25/Z+ . 0732422/Z3 
P1=1 .0+.1171875/Z2- .1441956/Z4 
Ol=.375/Z- .10253906/Z3 
X0»(Z-.25*PI) 

X1=(Z-.75*PI) 

E0=CEXP(-j*X0) 

E1=CEXP(-j*X1) 

AM'CSQRT { P2/Z ) 

H0=AM*(P0-j*OO)*E0 
Hl=AH*(P1 - j*01 )*E1 

END  IF 

RETURN 

END 


*******1 


r********************** 


SUBROUTINE  FACTOR  <A,P,N,NMX) 

PERFORMING  LU-OECOMPOSI TION  WITH  PIVOTING  ON  THE  A-MATRIX 

INTEGER  I , J, JP1,K,P(NMX),PJ,PR,R,RM1,RP1 
REAL  OMAX.ELMAG 
COMPLEX  A(NMX,NMX>. 0(120) 

DO  60  R«1,N 

DO  10  K«1,N 
10  OOC)=A((C,R) 

RM1*R-1 

IF(RMI.LT.I)  GO  TO  31 
DO  30  J»1 ,RM1 
PJ«P(J) 

A(J,R)«0(PJ> 

0(PJ)*0(J> 

JP1*J+1 
DO  20  I*JP1 ,N 


101 


20  D<l)*D<I)-A(I,J)*A(J,R> 

30  CONTINUE 

31  CONTINUE 
ONAX*CABS<D(R» 

P<R)*R 

RP1-R+1 

IF(RPI.GT.N)  GO  TO  41 
DO  40  I-RP1.N 

ELMAG*CABS(D(I)) 
lF(ELMAG.LT.DMAX)  GO  TO  40 
DNAX-ELMAG 

40  P<R)*I 

41  CONTINUE 
PR*P(R) 

A(R,R>«D(PR) 

D(PR)*D(R) 

IF(RPI.GT.N)  GO  TO  51 
00  50  I*RP1,N 

50  A<t,R)*0<I)/A<R,R> 

51  CONTINUE 
60  CONTINUE 

RETURN 

END 


C 

g  A****************- 

c 


SUBROUTINE  SOLVE  (A,B,P,N,NHX) 

C  BACK  SUBSTITUTION  TO  INVERT  THE  LINEAR  SYSTEM 
COMPLEX  A(NMX,NMX),B(NMX),Y(120),SUH 
INTEGER  I,IP1,J,K,P{NMX),PPI 
DO  20  1*1, N 

PPI*P(I) 

Y(I)*B<PPI) 

B<PPI)*B(I) 

IP1-I-H 

IF(IPI.GT.N)  GO  TO  11 
DO  10  4-IP1.N 

B( J)*B( J)-A( J , I >*Y( I ) 

10  CONTINUE 

11  CONTINUE 
20  CONTINUE 

DO  40  K*1,N 

I*N*K+1 

SUM*(0.,0.) 

IP1*H1 

IF(IPI.GT.N)  GO  TO  31 
DO  30  J*IP1 ,N 

SUM*SUM+AU,J)*B(J) 

30  CONTINUE 

31  CONTINUE 
B(I)*(Y(I)-SUM)/A(I,I) 

40  CONTINUE 
RETURN 
END 
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APPENDIX  E 


SAMPLE  INPUT/OUTPUT  DATA  FILES  USED  WITH  MONO 

This  Appendix  lists  the  input  and  output  data  file 
structures  used  within  MONO.  There  are  two  different  types 
of  input  data  files,  one  that  is  entered  at  the  keyboard  at 
run  time  and  one  which  allows  unattended  processing  of  a 
number  of  runs  without  user  intervention.  In  each  case,  the 
data  fields  may  be  either  arranged  in  column  format  with  one 
entry  per  line,  or  all  on  one  line  separated  by  commas. 
Additionally,  a  variable  type  is  set  by  the  program  which 
allows  a  real  variable  to  read  both  an  input  of  "1"  or  "1." 
as  the  same. 

A.  INPUT  DATA  FILE  STRUCTURE 

This  data  file  may  be  adjusted  to  allow  the  user  the 
option  of  selecting  either  the  output  file  name  or,  for  the 
case  of  M0N08 ,  the  calculation  of  the  electric  field  along  the 
antenna  at  run  time.  Below  is  a  commented  data  set  that  could 
be  used.  Note  the  comments  are  NOT  to  be  entered  in  the  data 
set  as  errors  would  occur. 
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Data  File 
Entry 

w 


299.8 


.24 


.01 

.02 

2.17 

30 

30 

60 

30 


d:\file.out 


Comment 

This  identifies  the  basis  for  all  physical 
measurements  of  the  antenna.  Using  a  "w" 
or  "W"  signifies  all  measurements  are  in 
terms  of  wavelength,  while  a  "m"  or  "M" 
indicates  measurements  are  in  terms  of 
meters . 

Driving  source  frequency  in  MHz 

Antenna  height,  k.  Can  be  entered  as 
"0 . 24” . 

Antenna  radius,  a. 

Gap  distance,  q. 

Ground  plate  spacing,  /. 

Number  of  modes  for  Region  I,  Nt. 

Number  of  modes  for  Region  II,  N2. 

Number  of  modes  for  Gap  expansion,  M^ 

Not  used  for  M0N08. 

Number  of  points  to  calculate  the  current 
along  the  antenna. 

DOS  path  and  name  to  assign  to  output 
file.  Maximum  length  is  25  characters. 
This  entry  may  also  be  entered  at  the  key¬ 
board  at  run  time. 


t»yM  or  i»nw  Indicates  whether  to  calculate  the  elec¬ 
tric  field  along  the  antenna.  Test  is  NOT 
case  sensitive.  If  yes,  data  is  stored  as 
' EZ . DAT '  in  the  current  directory  in  for¬ 
mat  for  PLOT. FOR  found  in  Appendix  B. 

This  may  be  entered  via  the  keyboard  at 
run  time. 


The  same  data  file  can  be  arranged  as  follows: 

W, 299.8, .24, . 01. 02, 2. 17, 30, 30, 60, 30, d:\file. out, n 
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The  data  file  is  easily  modified  to  allow  the  user  to 
group  a  series  of  runs  together  to  be  batch  processed  later 
using  the  standard  DOS  redirection  command.  For  example,  to 
use  input  data  file  named  'RUNl1  and  to  have  the  normal  screen 
output  redirected  to  a  running  log  file  named  'OUTPUT.LOG', 
the  correct  syntax  would  be: 

MONO  < [path ] RUNl  » [ path ] OUTPUT .  LOG 
where  path  is  the  DOS  path  where  the  individual  files  are 
stored.  This  parameter  is  optional  and  defaults  to  the 
current  directory.  A  number  of  runs  could  then  be  completed 
without  user  input  by  combining  a  series  of  these  commands  in 
one  batch  file.  The  file  structure  is  the  same  as  that  listed 
above  except  that  an  'n'  is  inserted  at  the  beginning  of  the 
file  to  indicate  that  a  separate  input  file  is  not  to  be  used. 
Therefore,  the  same  data  file  used  above  would  be  modified  as 
seen  below: 

n,w, 299 . 8 , .24, .01.02,2. 17, 3 0,30, 60, 30, d:\file. out, n 

B.  SAMPLE  OUTPUT  FILE 

The  output  file  generated  by  M0N07  for  the  inputs  used 
above  is  listed  below:  (*  indicates  not  included  in  output 
from  MONOS) 


; 
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Driving  Frequency, 

299.800  MHZ 

Monopole  Height, 

.240  meters 

(  .240 

Lambda) 

Monopole  Radius, 

.010  meters 

(  .010 

Lambda) 

Feed  Gap  Distance, 

.010  meters 

(  .010 

Lambda) 

Upper  Plane  Height, 

2.171  meters 

(2.170 

Lambda) 

Region  I  Coeffs.  (Nl) : 

30 

Region  II  Coeffs.  (N2) : 

30 

No.  Gap  Coeffs.  (Ml) : 

60 

No.  Points  for  I(z): 

30 

Listing  the  Expansion  Coefficients 
N  and  B(N) : 

0  (-9 . 259474E-05 ,  -3 . 712934E-04 ) 

1  (-1 . 716743E-04 , -7 . 285348E-04 ) 

2  (— 1. 250944E-04 ,  -6 . 843712E-04) 

3  (-8 . 932992E-06 , -5 . 970440E-04 ) 

4  (6 . 372479E-04 ,  -3 . 034107E-04) 

5  (-7 . 747178E-04 ,  -5 . 976821E-04) 

6  (-4. 176759E-04,-5. 106181E-04) 

7  (-2 . 877242E-04 ,  -4 . 189367E-04 ) 

8  (-2 . 029108E-04 ,  -3 . 268855E-04 ) 

9  (— 1. 371401E-04 ,  -2 . 384746E-04) 

10  (-8 . 359972E-05 , -1. 573042E-04 ) 

11  (-4 . 017067E-05, -8 . 621241E-05) 

12  (-6 . 131331E-06 , -2 . 729695E-05) 

13  (1 . 893293E-05 , 1 . 824955E-05) 

14  (3 . 552895E— 05, 5 . 013411E-05) 

15  (4 . 440288E-05 , 6 . 891159E-05) 

16  (4 . 656793E-05 , 7 . 587663E-05) 

17  (4 . 326494E-05 , 7 . 290936E-05) 

18  (3 . 588156E-05 , 6. 228465E-05) 

19  (2 . 58803 IE-05 , 4 . 649168E-05 ) 

20  ( 1 . 465865E-05 , 2 . 799412E-05) 


21  (3 . 497886E— 06 ,9. 087322E-06) 

22  (-6 . 533222E-06 , -8 . 273393E-06) 

23  (-1. 463069E-05 , -2 . 258123E-05) 

24  (-2 . 028524E-05 , -3 . 284282E-05) 

25  (-2 . 328759E-05, -3 . 859688E-05) 

26  (-2 . 370753E-05, -3 . 988246E-05) 

27  (-2 . 185403E-05, -3 .717663E-05) 

28  (-1 . 821469E-05 , -3 . 127903E-05) 

29  (-1. 339421E~05, -2 . 322528E-05) 

30  (-8 . 037723E-06, -1. 412969E-05) 
Listing  the  Current  on  the  Monopole 


Position 

Distance 

1(2)  (A) 

Number 

(Meters) 

Mag 

1 

.010 

2.316E-02 

-15.04 

2 

.018 

2 . 309E-02 

-15.98 

3 

.026 

2 . 299E-02 

-17.24 

4 

.034 

2 . 284E-02 

-18.58 

5 

.042 

2 . 262E-02 

-19.80 

6 

.050 

2.234E-02 

-20.78 
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7 

.058 

2.2f.  Jt,  '2 

-21.49 

8 

.066 

2 . 161E-02 

-21.97 

9 

.073 

2 . 121E-02 

-22.34 

10 

.081 

2 . 082E-02 

-22.72 

11 

.089 

2 . 044E-02 

-23.19 

12 

.097 

2 . 007E-02 

-23.78 

13 

.105 

1.970E-02 

-24.44 

14 

.113 

1.929E-02 

-25.11 

15 

.121 

1.881E-02 

-25.70 

16 

.129 

1. 822E-02 

-26.16 

17 

.137 

1.752E-02 

-26.50 

18 

.145 

1.669E-02 

-26.74 

19 

.153 

1. 574E-02 

-26.95 

20 

.161 

1.467E-02 

-27.19 

21 

.169 

1 . 349E-02 

-27.49 

22 

.177 

1 . 222E-02 

-27.86 

23 

.185 

1. 088E-02 

-28.27 

24 

.193 

9 . 485E-03 

-28.64 

25 

.200 

8 . 069E-03 

-28.91 

26 

.208 

6.671E-03 

-29.03 

27 

.216 

5.331E-03 

-28.99 

28 

.224 

4 . 091E-03 

-28.84 

29 

.232 

2.987E-03 

-28.72 

30 

.240 

2 . 042E-03 

-28.84 

Input  Resistance  = 

41.703650 

Input  Reactance  = 

11.206850 
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APPENDIX  F 


SAMPLE  INPUT/OUTPUT  DATA  FILES  USED  WITH  NEC 


This  Appendix  lists  an  input  data  file  used  NEC  and  the 
corresponding  output  file.  The  run  selected  was  for  an 
antenna  with  the  following  physical  characteristics  using  one 
reflection: 


Driving  Frequency 
Antenna  Height 
Antenna  Radius 
Gap  Distance 
Ground  Plate  Spacing 


299.8  MHz 
0.48  wavelengths 
0.01  wavelengths 
0.02  wavelenghts 
1.4  wavelengths 


CE  Run  11:  Honopole  Antenna  of  Height  .48, Radius  .01,. 02, 1.40 

GW  1,  24,  0,  0,  0,  0,  0,  .48,  .01 

GW  2,  48,  0,  0,  2.32,  0,  0,  3.28,  .01 

GE  1 

GN  1 

EK 

EX  0  1  1  0  1.0 

EX  0  2  24  0  -1.0 

EX  0  2  25  0  -1.0 

XQ 
EN 
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Sample  output  file  from  the  above  data  is  shown  below: 


r*******************************i 

NUMERICAL  ELECTROMAGNETICS  CODE 


- COMMENTS - 

Run  11:  Monopola  Antatma  of  Haight  .*8, Radius  .01,. 02,1. 10 


-  -  -  STRUCTURE  SPECIFICATION - 

COORDINATES  MUST  BE  INPUT  IN 
METERS  OR  BE  SCALED  TO  METERS 
BEFORE  STRUCTURE  INPUT  IS  ENDED 


WIRE  NO.  OF  FIRST 

LAST  TAG 


NO. 

SEG. 

XI 

NO. 

Yl 

Z1 

X2 

Y2 

Z2 

RADIUS 

SEG. 

SEG. 

1 

21 

. 00000 

1 

. 00000 

.00000 

.00000 

.00000 

.  18000 

.01000 

21 

1 

2 

72 

.00000 

2 

.00000 

2.32000 

.00000 

.00000 

3.28000 

.01000 

18 

25 

GROUND  PLANE  SPECIFIED. 

WHERE  WIRE  ENDS  TOUCH  GROUND,  CURRENT  WILL  BE  INTERPOLATED  TO  IMAGE  IN  GROUND  PLANE. 


TOTAL  SEGMENTS  USED-  72  NO.  SEG.  IN  A  SYfWETRIC  CELL-  72  SYMMETRY  FLAG-  0 


-  MULTIPLE  WIRE  JUNCTIONS  - 
JUNCTION  SEGMENTS  (-  FOR  END  1,  FOR  END  2) 
NONE 


-  SEGMENTATION  DATA  - 

COORDINATES  IN  METERS 

1+  AND  I-  INDICATE  THE  SEGMENTSBEFCRE  AND  AFTER  I 


SEG. 

COORDINATES  OF  SEG. 

CENTER 

SEG. 

ORIENTATION  ANGLES 

WIRE 

CONNECTION  DATA 

TAG 

NO. 

X 

Y 

z 

LENGTH 

ALPHA 

BETA 

RADIUS 

I- 

I 

1+ 

NO. 

1 

.00000 

.00000 

.01000 

.02000 

90.00000 

.00000 

.01000 

1 

1 

2 

1 

2 

.00000 

.00000 

.03000 

.02000 

90.00000 

. 00000 

.01000 

1 

2 

3 

I 

3 

.00000 

.00000 

.05000 

.02000 

90 . 00000 

.00000 

.01000 

2 

3 

4 

1 

1 

.00000 

.00000 

.07000 

.02000 

90.00000 

.00000 

.01000 

3 

4 

5 

1 

5 

.00000 

.00000 

.09000 

.02000 

90.00000 

.00000 

.01000 

* 

5 

6 

1 

6 

.00000 

.00000 

.11000 

.02000 

90.00000 

.00000 

.01000 

5 

6 

7 

1 

7 

.00000 

.00000 

. 13000 

.02000 

90.00000 

.00000 

.01000 

6 

7 

8 

1 

8 

.00000 

.00000 

. 15000 

.02000 

90.00000 

.00000 

.01000 

7 

8 

9 

1 

9 

.00000 

.00000 

.  17000 

.02000 

90.00000 

.00000 

.01000 

8 

9 

10 

1 

10 

.00000 

.00000 

.  19000 

.02000 

90.00000 

.00000 

.01000 

9 

10 

11 

1 

11 

.00000 

.00000 

.21000 

.02000 

90.00000 

.00000 

.01000 

10 

11 

12 

1 

12 

.00000 

.00000 

.23000 

.02000 

90.00000 

. 00000 

.01000 

11 

12 

13 

I 

13 

.00000 

.00000 

.25000 

.02000 

90.00000 

.00000 

.01000 

12 

13 

14 

1 

109 


14 

.00000 

.00000 

.27000 

.02000 

90.00000 

. 00000 

.01000 

13 

14 

15 

1 

IS 

.00000 

.00000 

.29000 

.02000 

90.00000 

.00000 

.01000 

14 

15 

16 

1 

16 

.00000 

.00000 

.31000 

. 02000 

90.00000 

.00000 

.01000 

15 

16 

17 

1 

17 

.00000 

.00000 

.33000 

.02000 

90 . 00000 

.00000 

.01000 

16 

17 

18 

1 

18 

.00000 

.00000 

.35000 

.02000 

90.00000 

.00000 

.01000 

17 

18 

19 

1 

19 

. 00000 

.00000 

.37000 

.02000 

90.00000 

.00000 

.01000 

18 

19 

20 

1 

20 

.00000 

.00000 

.39000 

.02000 

90.00000 

.00000 

.01000 

19 

20 

21 

1 

21 

.00000 

.00000 

.41000 

. 02000 

90.00000 

.00000 

.01000 

20 

21 

22 

1 

22 

.00000 

.00000 

.43000 

.02000 

90.00000 

.00000 

.01000 

21 

22 

23 

1 

23 

. 00000 

.00000 

.45000 

.02000 

90.00000 

.00000 

.01000 

22 

23 

24 

1 

24 

.00000 

.00000 

.47000 

.02000 

90.00000 

.00000 

.01000 

23 

24 

0 

1 

25 

.00000 

.00000 

2.33000 

.02000 

90.00000 

.00000 

.01000 

0 

25 

26 

2 

26 

.00000 

.00000 

2.35000 

.02000 

90.00000 

. 00000 

.01000 

25 

26 

27 

2 

27 

.00000 

.00000 

2.37000 

.02000 

90.00000 

.00000 

.01000 

26 

27 

28 

2 

26 

.00000 

.00000 

2.39000 

.02000 

90.00000 

.00000 

.01000 

27 

28 

29 

2 

29 

.00000 

.00000 

2.41000 

.02000 

90.00000 

.00000 

.01000 

28 

29 

30 

2 

30 

.00000 

.00000 

2.43000 

.02000 

90.00000 

.00000 

.01000 

29 

30 

31 

2 

31 

.00000 

.00000 

2.45000 

.02000 

90.00000 

. 00000 

.01000 

30 

31 

32 

2 

32 

.00000 

.00000 

2.47000 

.02000 

90.00000 

.00000 

.01000 

31 

32 

33 

2 

33 

.00000 

.00000 

2.49000 

.02000 

90.00000 

.00000 

.01000 

32 

33 

34 

2 

34 

.00000 

.00000 

2.51000 

.02000 

90.00000 

.00000 

.01000 

33 

34 

35 

2 

35 

.00000 

.00000 

2.53000 

.02000 

90.00000 

.00000 

.01000 

34 

35 

36 

2 

36 

.00000 

.00000 

2.55000 

.02000 

90.00000 

.00000 

.01000 

35 

36 

37 

2 

37 

.00000 

.00000 

2.57000 

.02000 

90.00000 

.00000 

.01000 

36 

37 

38 

2 

38 

.00000 

.00000 

2.59000 

.02000 

90.00000 

.00000 

.01000 

37 

38 

39 

2 

39 

.00000 

.00000 

2.61000 

.02000 

90.00000 

.00000 

.01000 

38 

39 

40 

2 

40 

.00000 

.00000 

2.63000 

.02000 

90.00000 

.00000 

.01000 

39 

40 

41 

2 

41 

.00000 

.00000 

2.65000 

.02000 

90.00000 

.00000 

.01000 

40 

41 

42 

2 

42 

.00000 

.00000 

2.67000 

.02000 

90.00000 

.00000 

.01000 

41 

42 

43 

2 

43 

.00000 

.00000 

2.69000 

.02000 

90 . 00000 

.00000 

.01000 

42 

43 

44 

2 

44 

.00000 

.00000 

2.71000 

.02000 

90.00000 

.00000 

.01000 

43 

44 

45 

2 

45 

.00000 

.00000 

2.73000 

.02000 

90.00000 

.00000 

.01000 

44 

45 

46 

2 

46 

.00000 

.00000 

2.75000 

.02000 

90.00000 

.00000 

.01000 

45 

46 

47 

2 

47 

.00000 

.00000 

2.77000 

.02000 

90.00000 

. 00000 

.01000 

46 

47 

48 

2 

46 

.00000 

.00000 

2.79000 

.02000 

90.00000 

.00000 

.01000 

47 

48 

49 

2 

49 

. 00000 

.00000 

2.81000 

.02000 

90.00000 

.00000 

.01000 

46 

49 

50 

2 

50 

.00000 

.00000 

2.83000 

.02000 

90.00000 

.00000 

.01000 

49 

50 

51 

2 

51 

.00000 

.00000 

2.85000 

.02000 

90.00000 

.00000 

.01000 

50 

51 

52 

2 

52 

. 00000 

.00000 

2.87000 

.02000 

90.00000 

.  00000 

.01000 

51 

52 

53 

2 

53 

.00000 

.00000 

2.89000 

.02000 

90.00000 

.00000 

.01000 

52 

53 

54 

2 

54 

.00000 

.00000 

2.91000 

.  02000 

90.00000 

.00000 

.01000 

53 

54 

55 

2 

55 

.00000 

.00000 

2.93000 

.02000 

90.00000 

.00000 

.01000 

54 

55 

56 

2 

56 

.00000 

.00000 

2.95000 

.02000 

90.00000 

.00000 

.01000 

55 

56 

57 

2 

57 

.00000 

. 00000 

2.97000 

.02000 

90.00000 

.00000 

.01000 

56 

57 

58 

2 

58 

.00000 

.00000 

2.99000 

.  02000 

90.00000 

.00000 

.01000 

57 

58 

59 

2 

59 

.00000 

.00000 

3.01000 

.02000 

90.00000 

.00000 

.01000 

58 

59 

60 

2 

60 

.00000 

.00000 

3.03000 

.02000 

90.00000 

.00000 

.01000 

59 

60 

61 

2 

61 

.00000 

.00000 

3.05000 

.02000 

90.00000 

.00000 

.01000 

60 

61 

62 

2 

62 

.00000 

.00000 

3.07000 

.02000 

90.00000 

.00000 

.01000 

61 

62 

63 

2 

63 

.00000 

.00000 

3.09000 

.  02000 

90.00000 

.00000 

.01000 

62 

63 

64 

2 

64 

.00000 

.00000 

3.11000 

.02000 

90.00000 

.  00000 

.01000 

63 

64 

65 

2 

85 

. 00000 

.00000 

3 . 13000 

.02000 

90.00000 

.00000 

.01000 

64 

65 

66 

2 

66 

.00000 

.00000 

3.15000 

.02000 

90.00000 

.00000 

.01000 

65 

66 

67 

2 

67 

.00000 

.00000 

3.17000 

.02000 

90.00000 

.  00000 

.01000 

66 

67 

68 

2 

68 

.00000 

.00000 

3 . 19000 

.02000 

90 . 00000 

.  00000 

.01000 

67 

68 

69 

2 

69 

.00000 

.00000 

3.21000 

.02000 

90.00000 

.  00000 

.01000 

68 

69 

70 

2 

70 

.00000 

.00000 

3.23000 

.02000 

90.00000 

.00000 

.01000 

69 

70 

71 

2 

71 

.00000 

.00000 

3.25000 

.02000 

90 . 00000 

.  00000 

.01000 

70 

71 

72 

2 

72 

.00000 

.00000 

3.27000 

.02000 

90.00000 

.  00000 

.01000 

71 

72 

0 

2 

O.OOOOOE+OO  0 .OOOOOE+OO 

GR 

1 

0 

0 

0  O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

*****  DATA  CARD  RO.  2 
O.OOOOOE+OO  0.00000E+00 

EK 

0 

0 

0 

0  0 . OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

*****  DATA  CARD  RO.  3 

0 . OOOOOE+OO  0 . 00000E+00 

EX 

0 

1 

1 

0  1. OOOOOE+OO 

O.OOOOOE+OO 

0 . OOOOOE+OO 

O.OOOOOE+OO 

*****  DATA  CARD  RO.  4 
O.OOOOOE+OO  O.OOOOOE+OO 

EX 

0 

2 

24 

0  -1. OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

*****  DATA  CARD  RO.  5 
O.OOOOOE+OO  O.OOOOOE+OO 

EX 

0 

2 

25 

0  -1. OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

110 


*****  DATA  CARO  NO.  6  XQ  0  0  0  0  0.00000E+00  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO 

0.00000E+00  O.OOOOOE+OO 


• - FREQUENCY - - 

FREQUENCY-  2.9980E+02  m2 
WAVELENGTH-  1.0000E+00  METERS 


APPROXIMATE  INTEGRATION  EMPLOYED  FOR  SEGMENTS  MORE  THAN  1.000  WAVELENGTHS  APART 
THE  EXTENDED  THIN  WIRE  KERNEL  WILL  BE  USED 


-  STRUCTURE  IMPEDANCE  LOADING  -  - 

THIS  STRUCTURE  IS  NOT  LOADED 


-  ANTENNA  ENVIRONMENT  - 

PERFECT  GROUND 


-  MATRIX  TIMING  - 

FILL-  1.500  MIN.,  FACTOR-  .36*  MIN. 


4 

-  ANTENNA  INPUT 

PARAMETERS  - 

-  - 

TAG 

(MHOS) 

SEG.  VOLTAGE 

POWER 

(VOLTS) 

CURRENT  (AMPS) 

IMPEDANCE 

(OHMS) 

ADMITTANCE 

NO. 

IMAG. 

NO.  REAL 

(WATTS) 

IMAG. 

REAL  IMAG. 

REAL 

IMAG. 

REAL 

1  1  1.00000E+00  O.OOOOOE+OO  2.*305*E-03  2.78753E-03  1.77699E+02-2.03799E+02  Z.*305*E~03 

2.7B753E-03  1.21527E-03 


2  48-1.00000E+00  O.OOOOOE+O0-2.*2866E-03-2.791*3E-03  1.77398E+02-2.03896E+02  2.42866E-03 

2.79143E-03  1.21A33E-03 

2  *9-1 . OOOOOE+OO  0 . OOOOOE+OO-2 . *2265E-03~2 . 79230E-03  1. 7727*E+02-2. 0*322E+02  2.*2265E-03 

2.79230E-03  1.21132E-03 


-  CURRENTS  AND  LOCATION  - 

DISTANCES  IN  WAVELENGTHS 


SEG. 

TAG 

COORD. 

OF  SEG. 

CENTER 

SEG. 

-  -  -  CURRENT  (AMPS)  - 

NO. 

NO. 

X 

Y 

z 

LENGTH 

REAL  IMAG. 

MAG. 

PHASE 

1 

1 

.0000 

.0000 

.0100 

.02000 

2.A305E-03  2.7875E-03 

3.698*E-03 

*8.91* 

2 

1 

.0000 

.0000 

.0300 

.02000 

2.*l*7E-03  1.0516E-03 

2.6338E-03 

23.53* 

3 

1 

.0000 

.0000 

.0500 

. 02000 

2. 3833E-03  -8.5797E-05 

2.38*8E-03 

-2.062 

* 

1 

.0000 

.0000 

.0700 

.02000 

2.3366E-03  -9.8162E-0* 

2.534*E-03 

-22.787 

3 

1 

.0000 

.0000 

.0900 

.02000 

2.275*E-03  -1.7685E-03 

2.8819E-03 

-37.856 

6 

1 

.0000 

.0000 

.1100 

.02000 

2.200*E-03  -2. *668E-03 

3.3055E-03 

-*8.267 

7 

1 

.0000 

.0000 

.1300 

.02000 

2. 1125E-03  -3.08*7E-03 

3. 7387E-03 

-55.595 

8 

1 

.0000 

.0000 

.1500 

.02000 

2. 0130E-03  -3.62*0E-03 

*.l*56E-03 

-60.950 

9 

1 

.0000 

.0000 

.1700 

.02000 

1 . 9031E-03  -*.08*lE-03 

* . S058E-03 

-65.016 

10 

1 

.0000 

.0000 

.1900 

.02000 

1.78*2E-03  -*.*630E-03 

* . 806SE-03 

-68.209 

11 

1 

.0000 

.0000 

.2100 

.02000 

1 . 6S79E-03  -*. 7589E-03 

5.0394E-03 

-70.792 

12 

1 

.0000 

.0000 

.2300 

.02000 

1.5258E-03  -*.9697E-03 

5. 1987E-03 

-72.932 

13 

1 

.0000 

.0000 

.2500 

.02000 

1.389SE-03  -5.09*3E-03 

5.2806E-03 

-7*.  7** 

1* 

1 

.0000 

.0000 

.2700 

.02000 

1.2508E-03  -5. 1325E-03 

5, 2827E-03 

-76.30* 

15 

1 

.0000 

.0000 

.2900 

.02000 

1.1113E-03  -5.08*lE-03 

3. 20* IE-03 

-77.670 

16 

1 

.0000 

.0000 

.3100 

.02000 

9. 7278E-0*  -*.9502E-03 

5.0**9E-03 

-78.882 

111 


17 

1 

.0000 

.0000 

.3300 

.02000 

8 . 3692E-04 

*4 . 7327E-03 

4 . 8061E-03 

-79.872 

18 

1 

.0000 

.0000 

.3500 

.02000 

7 . 0528E-04 

-4.4340E-03 

4.4897E-03 

-80.962 

19 

1 

.0000 

.0000 

.3700 

. 02000 

5.7941E-04 

-4.0569E-03 

4 . 0981E-03 

-81.872 

20 

1 

.0000 

.0000 

.3900 

.02000 

4.6064E-04 

-3.6044E-03 

3 . 6337E-03 

-82.717 

21 

1 

.0000 

.0000 

.4100 

.02000 

3.S016E-04 

-3.0785E-03 

3.0983E-03 

-83.510 

22 

1 

.0000 

.0000 

.4300 

.02000 

2.4894E-04 

-2.4785E-03 

2.4909E-03 

-84.264 

23 

1 

.0000 

.0000 

.4500 

.02000 

1.5733E-04 

-1.7963E-03 

1.8032E-03 

-84.995 

24 

1 

.0000 

.0000 

.4700 

.02000 

7.2003E-05 

-9.5821E-04 

9.6091E-04 

-85.703 

25 

2 

.0000 

.0000 

2.3300 

.02000 

-8.2718E-05 

9.S969E-04 

9. 6325E-04 

94 . 926 

26 

2 

.0000 

.0000 

2.3500 

.02000 

-1.7730E-04 

1.7990E-03 

1.8077E-03 

95.629 

27 

2 

.0000 

.0000 

2.3700 

.02000 

-2.7635E-04 

2.4819E-03 

2 . 4972E-03 

96.353 

28 

2 

.0000 

.0000 

2.3900 

.02000 

-3.8402E-04 

3 . 0826E-03 

3 . 1064E-03 

97.101 

29 

2 

.0000 

.0000 

2.4100 

.02000 

-5.0002E-04 

3.6090E-03 

3.6434E-03 

97.888 

30 

2 

.0000 

.0000 

2.4300 

.02000 

-6.2347E-04 

4.0618E-03 

4 . 1094E-03 

98.727 

31 

2 

.0000 

.0000 

2.4500 

.02000 

-7.5316E-04 

4.4391E-03 

4 . 5025E-03 

99.629 

32 

2 

.0000 

.0000 

2.4700 

.02000 

-8.8769E-04 

4 . 7379E-03 

4.8203E-03 

100.612 

33 

2 

.0000 

.0000 

2.4900 

.02000 

-1.0256E-03 

4 . 9554E-03 

5.0804E-03 

101.693 

34 

2 

.0000 

.0000 

2.5100 

.02000 

-1.1652E-03 

5.0891E-03 

S.2208E-03 

102.896 

35 

2 

.0000 

.0000 

2.5300 

.02000 

-1.3049E-03 

5. 1373E-03 

5.3004E-03 

104 . 252 

36 

2 

.0000 

.0000 

2.5500 

.02000 

-1.4430E-03 

S.0989E-03 

5.2992E-03 

105.601 

37 

2 

.0000 

.0000 

2.5700 

. 02000 

-1.5778E-03 

4 . 9737E-03 

S.2180E-03 

107.600 

38 

2 

.0000 

.0000 

2.5900 

.02000 

-1.7076E-03 

4 . 7623E-03 

5.0592E-03 

109.726 

39 

2 

.0000 

.0000 

2.6100 

.02000 

-1.8308E-03 

4 . 4659E-03 

4.8266E-03 

112.292 

40 

2 

.0000 

.0000 

2.6300 

.02000 

-1.9460E-03 

4.0864E-03 

4 . 5261E-03 

115.464 

41 

2 

.0000 

.0000 

2.6500 

.02000 

-2.0515E-03 

3.6256E-03 

4 . 1658E-03 

119.502 

42 

2 

.0000 

.0000 

2.6700 

.02000 

-2.1461E-03 

3.08S6E-03 

3.7585E-03 

124 . 820 

43 

2 

.0000 

.0000 

2.6900 

.02000 

-2.2286E-03 

2.4669E-03 

3.3245E-03 

132.094 

44 

2 

.0000 

.0000 

2.7100 

.02000 

-2.2978E-03 

1.7678E-03 

2 . 8983E-03 

142.427 

45 

2 

.0000 

.0000 

2.7300 

.02000 

-2.3532E-03 

9.8018E-04 

2.5492E-03 

157.387 

46 

2 

.0000 

.0000 

2.7500 

.02000 

-2.3937E-03 

8 . 3523E-05 

2.39S2E-03 

178.002 

47 

2 

.0000 

.0000 

2.7700 

.02000 

-2.4190E-03 

-1.0547E-03 

2.6389E-03 

-156.441 

48 

2 

.0000 

.0000 

2.7900 

.02000 

-2.4287E-03 

-2.7914E-03 

3.7001E-03 

-131.025 

49 

2 

.0000 

.0000 

2.8100 

.02000 

-2.4226E-03 

-2.7923E-03 

3.6968E-03 

-130.945 

50 

2 

.0000 

.0000 

2.8300 

.02000 

-2.4010E-03 

-1.0573E-03 

2.6235E-03 

-156.235 

51 

2 

.0000 

.0000 

2.8500 

.02000 

-2.3641E-03 

7.9330E-05 

2.3655E-03 

178.078 

52 

2 

.0000 

.0000 

2.8700 

. 02000 

-2.3124E-03 

9.7436E-04 

2.5093E-03 

157.151 

53 

2 

.0000 

.0000 

2.8900 

.02000 

-2.2466E-03 

1.7605E-03 

2.8542E-03 

141.917 

54 

2 

.0000 

.0000 

2.9100 

.02000 

-2.1675E-03 

2.4580E-03 

3.2772E-03 

131.406 

55 

2 

.0000 

.0000 

2.9300 

.02000 

-2.0762E-03 

3.0752E-03 

3.7104E-03 

124.025 

56 

2 

.0000 

.0000 

2.9500 

.02000 

-1.9738E-03 

3.6139E-03 

4.117BE-03 

118.642 

57 

2 

.0000 

.0000 

2.9700 

.02000 

-1.8817E-03 

4.0734E-03 

4.4786E-03 

114.562 

58 

2 

.0000 

.0000 

2.9900 

.02000 

-1.7413E-03 

4.4518E-03 

4.7803E-03 

111.362 

59 

2 

.0000 

.0000 

3.0100 

.02000 

-1.6141E-03 

4.7472E-03 

5.0141E-03 

108.779 

60 

2 

.0000 

.0000 

3.0300 

.02000 

'1 . 4819E-03 

4.9578E-03 

5.1745E-03 

106.641 

61 

2 

.0000 

.0000 

3.0500 

.02000 

-1.3461E-03 

5.0823E-03 

5.2575E-03 

104.834 

62 

2 

.0000 

.0000 

3.0700 

.02000 

-1.2084E-03 

5.1202E-03 

5.2609E-03 

103.280 

63 

2 

.0000 

.0000 

3.0900 

.02000 

-1.0707E-03 

5.0718E-03 

5.1836E-03 

101.921 

64 

2 

.0000 

.0000 

3.1100 

.02000 

-9.3452E-04 

4.9381E-03 

5.0257E-03 

100.716 

65 

2 

.0000 

.0000 

3.1300 

.02000 

-8.0147E-04 

4.7209E-03 

4.7884E-03 

99.635 

66 

2 

.0000 

.0000 

3.1500 

.02000 

-6.7310E-04 

4.4227E-03 

4.4736E-03 

98.654 

67 

2 

.0000 

.0000 

3.1700 

.02000 

-5.5087E-04 

4.0463E-03 

4.0836E-03 

97.753 

68 

2 

.0000 

.0000 

3 . 1900 

.02000 

-4.3609E-04 

3 . 5946E-03 

3.6210E-03 

96.917 

69 

2 

.0000 

.0000 

3.2100 

.02000 

-3.2990E-04 

3 . 0698E-03 

3.0875E-03 

96.134 

70 

2 

.0000 

.0000 

3.2300 

.02000 

-2.3317E-04 

2.4712E-03 

2.4822E-03 

95.390 

71 

2 

.0000 

.0000 

3.2500 

.02900 

-1.4630E-04 

1. 7909E-03 

1.7969E-03 

94.670 

72 

2 

.0000 

.0000 

3.2700 

.02000 

-6.6339E-05 

9.5517E-04 

8.5747E-04 

93.973 

-  TOWER  BUDGET  - 

INPUT  TOWER  -  3. 6409E-03  WATTS 
RADIATED  POWER-  3.6409E-03  WATTS 
STRUCTURE  LOSS-  O.OOOOE+OO  WATTS 
NETWORK  LOSS  -  O.OOOOE+OO  WATTS 
EFFICIENCY  -  X00.00  PERCENT 


*****  DATA  CARD  NO.  7  EN  0  0  0  0  O.OOOOOE+OO  0.00000E+00  O.OOOOOE+OO  O.OOOOOE+OO 

O.OOOOOE+OO  O.OOOOOE+OO 

RUN  TIME  -  1.952 
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